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Abstract 

An analytical, macroscopic approach to SU(N) Yang-Mills thermodynam- 
ics is developed. This approach self-consistently assumes that at a tempera- 
ture much larger than the Yang-Mills scale Aym,n (embedded and noninteract- 
ing) SU(2) calorons of trivial holonomy form an adjoint Higgs field (electric 
phase). Macroscopically, this field turns out to be thermodynamically and 



quantum mechanically stabilized. As a consequence, the problem of the in- 
frared instability in the perturbative loop expansion of thermodynamical po- 
tentials, generated by the soft magnetic modes, is resolved. An evolution equa- 
tion with two fixed points follows for the effective gauge coupling e{T) from 
self-consistent thermodynamics involving the ground-state and its quasiparti- 
cle excitations. A plateau value of e(T), which is an attractor of the evolution, 
is consistent with the existence of isolated magnetic monopoles of conserved 
charge being generated by dissociating calorons of nontrivial holonomy. The 
(up to negligible corrections exact) one- loop and downward evolution of e{T) 
predicts the condensation of magnetic monopoles in a 2"''-order phase tran- 
sition at a critical temperature Te,c- At Te,c tree-level massive gauge modes 
decouple thermodynamically. This is the confinement phase transition identi- 
fied in lattice simulations. For N=2 we compute the critical exponent taking 
the mass of the dual photon as an order parameter. For arbitrary N we show 
the restoration of the global electric symmetry in the monopole condensed 
(magnetic) phase by investigating the Polyakov loop in the effective theory. 
The magnetic gauge coupling g{T) starts its downward evolution from zero at 
Te,c and runs into a logarithmic pole at Tm,c < Te,c- At Tm,c center-vortex 
loops condense, the abelian gauge modes decouple thermodynamically, and 
the equation of state is p = —P (zero entropy density). The Hagedorn tran- 
sition to the vortex condensing phase (center phase) goes with a complete 
breakdown of the local magnetic Z^j symmetry. After a rapid reheating in 
terms of (intersecting) center-vortex loops has taken place the ground-state 
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pressure vanishes identically on tree level. This result is protected against 
radiative corrections. Throughout the electric and the magnetic phase and 
for N=2,3 we compute the temperature evolution of the (infrared sensitive) 
pressure and energy density and for the (infrared insensitive) entropy density 
and compare our results with lattice data. We show that the disagreement for 
the two former quantities at low temperature (negative pressure) originates 
from severe finite-size artefacts in lattice simulations. For the entropy density 
we obtain excellent agreement with lattice results. The implications of our 
results for particle physics and cosmology are discussed. 
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1 Introduction 



The beauty and usefulness of the gauge principle for local field theories is generally 
appreciated. Yet, in a perturbative approach to gauge theories like the SM and 
its (non)supersymmetric extensions it is hard if not impossible to convincingly ad- 
dress a number of recent experimental and observational results in particle physics 
and cosmology: nondetection of the Higgs particle at LEP [T], indications for a 
rapid thermalization and strong collective behavior in the early stage of an ultra- 
relativistic heavy-ion collision |21 E], dark energy and dark matter at present, a 
strongly favored epoch of cosmological inflation in the early Universe jH El El El , 
and the likely existence of intergalactic magnetic fields [3 El- An analytical and 
nonperturbative approach to strongly interacting gauge theories may further our 
understanding of these phenomena. 

It is difficult to gain insights in the dynamics of a strongly interacting field 
theory by analytical means. We conjecture with Ref. ^| that a thermodynamical 
approach is an appropriate starting point for such an endeavor. On the one hand, 
this conjecture is reasonable since a strongly interacting system being in equilib- 
rium communicates perturbations almost instantaneously due to rigid correlations. 
Thus equilibrium is restored very rapidly. On the other hand, the requirement of 
thermalization poses strong constraints on the construction of a macroscopic, effec- 
tive theory for the ground state and its (quasiparticle) excitations. The objective 
of the present paper is the thermodynamics of SU(N) Yang-Mills theories in four 
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dimensions. 

Let us very briefly recall some aspects of the analytical approaches to thermal 
SU(N) Yang-Mills theory as they are discussed in the literature. Because of asymp- 
totic freedom ^] one would naively expect thermal perturbation theory to 
work well for temperatures T much larger than the Yang-Mills scale Aym,n since the 
gauge coupling constant g(T) logarithmically approaches zero for ^ — >■ oo. It 
is known for a long time that this expectation is too optimistic since at any tem- 
perature perturbation theory is plagued by instabilities arising from the infrared 
sector (weakly screened, soft magnetic modes P^J). As a consequence, the pressure 
P can be computed perturbatively only up to (and including) order g^. The effects 
of resummations of one- loop diagrams (hard thermal loops), which rely on a scale 
separation in terms of the small value of the coupling constant g, are summarized 
in terms of a nonlocal effective theory for soft and semi-hard modes ^1]. In the 
computation of radiative effects based on this effective theory infrared effects due to 
soft modes still appear in an uncontrolled manner. This has lead to the construction 
of an effective theory where soft modes are collectively described in terms of classical 
fields whose dynamics is influenced by integrated semi-hard and hard modes [1^1 CHI • 
In Quantum Chromodynamics (QCD) a perturbative calculation of P was pushed 
up to order g^log g, and an additive 'nonperturbative' term at this order was fitted 
to lattice results ^7]. Within the perturbative orders a poor convergence of the 
expansion is observed for temperatures not much larger than the MS scale. While 
the work in jT^ is a computational masterpiece it could, by definition, not shed 
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light on the nonperturbative physics of the infrared sector. Screened perturbation 
theory, which rehes on a separation of the tree-level Yang-Mills action using varia- 
tional parameters, is a very interesting idea. Unfortunately, this approach generates 
temperature dependent ultraviolet divergences in its presently used form, see [TH] 
for a recent review. 

The purpose of this paper is to report in a detailed way^ on a nonperturbative 
and analytical approach to the thermodynamics of SU(N) Yang-Mills theory (see 
PU] for intermediate stages). Conceptually, this approach is similar to the macro- 
scopic Landau-Ginzburg-Abrikosov (LGA) theory for superconductivity in metals 
\21\ 122] ■ Recall, that this theory does not derive the condensation of Cooper pairs 
from first principles but rather describes the condensate by a nonvanishing am- 
plitude of a complex scalar field (local order parameter) which is charged under 
the electromagnetic gauge group U(l). This nonvanishing amplitude is driven by 
a phenomenologically introduced potential V. As a consequence, a (macroscopic) 
U(l) gauge field Op, which is deprived of the microscopic gauge-field fluctuations 
associated with the formation of Cooper pairs and their subsequent condensation, 
acquires mass, the U(l) symmetry is spontaneously broken, and physical phenomena 
originating from this breakdown can be explored in dependence of the parameters 
appearing in the effective action, and in dependence of an external magnetic field 

and/or temperature. 

^Some aspects of the low-temperature physics are revised in the present paper as compared to 
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When applying this idea to the construction of a macroscopic theory for SU(N) 
Yang-Mills thermodynamics (YMTD) one is in a much better position as far as the 
uniqueness of the stabilizing potentials in each phase of the theory is concerned. 
These potentials are determined by thermodynamics and the requirement that, in a 
first step of the construction, they admit energy- and pressure-free macroscopic con- 
figurations describing the collective effects in an ensemble of energy- and pressure- 
free, noninteracting, and self-dual topological field configurations in the underlying 
theory. If a particular phase supports propagating gauge modes then, in a second 
step, the interactions between these topological defects are treated by solving the 
macroscopic gauge-field equations in terms of a pure gauge configuration in the 
background of the (inert) energy- and pressure-free scalar field. 

More specifically, we assume that at a large temperature a macroscopic adjoint 
scalar field (f) is generated by a dilute gas of trivial- holonomy calorons ^ fl^ . Calorons 
are Bogomolnyi-Prasad-Sommerfield (BPS) saturated (or self-dual) solutions |i25| to 
the classical Yang-Mills equations of motion in four-dimensional Euclidean space- 
time^ (time coordinate r is compactified on a circle, < r < t]^) with varying 
topological charge and embedding in SU(N). Calorons are topologically nontrivial, 
saturate the lowest possible value of the Euchdean action in a given topological 

sector, and thus are energy- and pressure-free configurations. Calorons with non- 
^We discuss in Sec. 12.51 why the critical temperature Tp for the onset of the formation of 

should be comparable to the cutoff-scale for the local field theory in four dimensions. 

■^Whenever we speak of a topological soliton this automatically includes the antisoliton. 
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trivial holonomy have BPS magnetic monopole constituents [221 I2H|- Their 
one-loop effective action scales with the three volume of the system jSO], and thus 
they should play no role in the thermodynamic limit. This conclusion, however, is 
no longer valid if the system generates domains of large but finite volume whose 
boundaries are generated by discontinuous changes of the color orientation of the 
field 0. Microscopically, nontrivial-holonomy calorons can be dynamically generated 
out of trivial-holonomy calorons by macroscopic domain collisions. These calorons 
dissociate into their magnetic monopole constituents subsequently, see for a 
discussion of the destabilizing effects of quantum fluctuations in the case of nontriv- 
ial holonomy. We thus anticipate the occurence of isolated magnetic charge whose 
abundance is governed by the (T dependent) typical volume of a domain. 

The property of vanishing energy and pressure of a caloron derives from its 
self-duality, that is, the kinetic and the interaction part in the Euclidean energy- 
momentum tensor precisely cancel when evaluated on a caloron. A potential Ve (the 
subscript E stands for electric phase) is constructed which stabilizes the modulus |0| 
for given T quantum mechanically and thermodynamically and which reflects the 
assumption that (p is composed of noninteracting, trivial-holonomy calorons. We 
would like to stress at this point that the effects of calorons are reflected as a 
dependence of the modulus of (p. As a consequence, the nontrivial-topology sector 
of the theory, indeed, is irrelevant at asymptotically large temperatures. 

A unique decomposition of each gauge-field configuration Ap contributing to the 
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partition function of the fundamental Yang-Mills theory is 

Ap = A^;^ + ap. (1) 

In Eq. dH) y4^°p is a minimally (that is, BPS saturated) topological part, represented 
by calorons, and ap denotes a remainder which has trivial topology. The configura- 
tions in having trivial holonomy would build up the ground state described by 
if no holonomy- changing interactions between them were allowed for. A change in 
holonomy by interactions, mediated by the topologically trivial sector, will macro- 
scopically manifest itself in terms of a finite, pure-gauge background a!'^ . A fluc- 
tuation 5ap about this background acquires mass by the adjoint Higgs mechanism 
if [0, 5ap] 7^ and thus the underlying gauge symmetry SU(N) is spontaneously 
broken to U(l)'^'^ at most. The degree of gauge symmetry breaking by calorons 
is a boundary condition set at an asymptotically high temperature Tp where the 
effect of oc T^^/^ on the Yang- Mills spectrum and its pressure ~ is very small 
since the ground state pressure scales as oc T . On the one hand, Higgs-mechanism 
induced masses provide infrared cutoffs in the loop expansions of thermodynami- 
cal quantities which resolves the problem of the infrared instability entcountered 
in perturbation theory. On the other hand, the compositeness scale |0| constrains 
the hardness of quantum fluctuations, and so the usual renormalization program 
needed to address ultra-violet divergences in perturbation theory is superfluous in 
the effective theory. Notice that this way of introducing a composite field (f) in an ef- 
fective description differs from the usual implementation of a Wilsonian flow, where 
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high-momentum modes are successively integrated out f T^ , since is built of 
calorons with an 'instanton' radius p not being smaller than \4>\'^- At the present 
stage the description of the ground state of an SU(N) Yang-Mills theory at high 
temperatures in terms of the field is self-consistent. The phase and the modulus 
of the field are derived from a microscopic definition in ^Ij . 

The nonperturbative approach to SU(N) YMTD proposed here implies the ex- 
istence of three rather than two phases: an electric phase at high temperatures, a 
magnetic phase for a small range of temperatures comparable to the scale Aym.n? 
and a center phase for low temperatures. The ground state in the magnetic phase 
confines fundamental, static test charges but allows for the propagation of massive, 
dual gauge bosons. The center phase is thermodynamically disconnected from the 
magnetic and the electric phase. In the electric phase an evolution equation for the 
effective gauge coupling constant e, which follows from the requirement of thermo- 
dynamical self-consistency of the one-loop expression for the pressure, has two fixed 
points associated with a highest and a lowest attainable temperature Tp and Te^c- 
It turns out that practically all strong-interaction effects of the theory are described 
by a temperature dependent ground-state pressure and tree-level masses for thermal 
quasiparticles such that higher loop corrections to thermodynamical quantities are 
tiny. 

At Tg; c the effective coupling e(T) exhibits a thin divergence of the form 

e(T)~-log(T-T,), (2) 
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and the theory undergoes a 2°'* order hke phase transition to a magnetic phase 
which is driven by the condensation of some of the magnetic monopoles residing 
inside dissociating nontrivial-holonomy calorons. In this transition a part of the 
continuous gauge symmetry, which survived the formation of the adjoint Higgs field 
in the electric phase, is broken spontaneously and the tree-level massive gauge 
modes of the electric phase decouple thermodynamically. In the case of submaximal 
gauge-symmetry breaking by in the electric phase condensates of magnetic and 
color-magnetic monopoles occur in the magnetic phase. The former are described 
by complex scalar and the latter by adjoint Higgs fields. In the case of maximal 
gauge symmetry breaking to U(l)'^"\ which we will only investigate in this paper, 
the (local) magnetic center symmetry .^N/2,mag and the continuous gauge symmetry 
U(l)'^/^"^ survive the transition to the magnetic phase, the (global) electric center 
symmetry .^N/2,eiec is fully restored. An approach to the thermodynamics in the 
magnetic phase, which is conceptually analoguous to the one in the electric phase, 
yields an evolution equation for the magnetic gauge coupling g which has two fixed 
points at Te,c and Tm,c (highest and lowest attainable temperature). Approaching 
Tm,c from above, the equation of state is increasingly dominated by the ground state 
contributions. At Tm,c we have 

p^-P. (3) 

The theory undergoes a phase transition to a phase whose ground-state is a con- 
densate of center-vortex loops. In this phase .^N.mag is entirely broken, and all gauge 
boson excitations are thermodynamically decoupled. Once each of the vortex-loop 
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condensates, described by nonlocally defined complex scalar fields, has relaxed to 
the one of the N degenerate minima of its potential, the energy density and the 
pressure of the ground state are precisely zero (no radiative corrections), and the 
system has created particles by local phase shifts of each vortex-condensate field 
which are associated with localized (intersecting) center fluxes forming closed loops. 
The corresponding density of states is over-exponentially rising implying that the 
magnetic-center transiton is of the Hagedorn type and thus nonthermal. 

There are many claims in the scenario outlined above. We will, step by step, 
verify them as we proceed. The paper is organized as follows: 

In Sec. 12] we explain our approach to the electric phase. We start with the basic 
assumption that it is noninteracting trivial-holonomy calorons that form a macro- 
scopic adjoint Higgs field (j) high temperatures (electric phase) and explore its 
consequences. A nonlocal definition for (j) is given. We then elucidate the details of 
the ground-state dynamics and the properties of topology-free gauge modes. Subse- 
quently, an evolution equation for the effective gauge coupling constant e is derived 
and solved, interpretations of the solution are given, and an argument is provided 
why the temperature Tp for the onset of caloron 'condensation' has to be compa- 
rable to the cutoff-scale for the local field-theory description in four dimensions. 
In a next step, we perform the counting of isolated magnetic monopoles species in 
the effective theory for the electric phase when assuming maximal gauge symme- 
try breaking by 0. The next part of SecE] is devoted to a discussion of two-loop 
corrections to thermodynamical quantities. For the SU(2) case we investigate the 
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simplest one-loop contribution to the 'photon' polarization and perform formal weak 
and strong coupling limits of this expression. We also discuss the implementation of 
thermodynamical self-consistency when higher loop corrections to the pressure are 
taken into account. 

In SecE] we investigate the magnetic phase, again assuming maximal gauge 
symmetry by caloron 'condensation': the pattern of gauge symmetry breaking by 
monopole condensation is explored, the thermodynamics of the ground state and 
its excitations is elucidated, an evolution equation for the magnetic gauge coupling 
constant g is derived. Solutions to this equations are obtained numerically and their 
implications are discussed. Finally, we discuss the Polyakov loop in the electric and 
the magnetic phase and compute the critical exponent of the phase transition for 
SU(2). 

In Sec.|3] we investigate the center phase. A nonlocal definition for the local 
fields describing the condensed center-vortex loops is given, their transformation 
properties under magnetic center rotations are determined, and their dynamics is 
discussed. 

In Sec. we derive a matching condition for the mass scales and A^ which 
appear in the respective potentials for the caloron and magnetic monopole conden- 
sates. 

In Sec. El we compute the temperature evolution of the thermodynamical poten- 
tials pressure, energy density, and entropy density throughout the electric and the 
magnetic phases at one loop for N=2,3 and compare our results with lattice data. 
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A conclusion and a discussion of likely implications of our results for particle 
physics and cosmology are given in Sec.[7| 

2 The electric phase 
2.1 Conceptual framework 

Our analysis is based on the following assumption about the ground-state physics 
characterizing SU(N) YMTD at high temperatures. 

At a temperature Tp ^ Aym,n SU(N) YMTD, defined on a Euclidean, four-dimensional, 
and fiat spacetime, generates an adjoint Higgs field out of noninteracting (dilute), 
trivial-holonomy SU(2) calorons. 

Calorons are BPS saturated solutions to the Euclidean equations of motion of SU(N) 
Yang- Mills theory at finite temperature |21l2ni IZ3 123 • One distinguishes SU(2) 
calorons according to their holonomy, that is, the behavior of the Polyakov loop 



at |f I — > cxD when evaluated on the solution. In Eq. (jH) V denotes the path-ordering 
symbol and g the gauge coupling constant of the SU(2) Yang- Mills theory. Trivial 
(nontrivial) holonomy means that have P|:j|^oo = 1 (7^ !)• In the former case the 
SU(2) caloron has no isolated magnetic-monopole constituents, in the latter case 
it exhibits a monopole and its antimonopole. The masses of these constituents are 
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Figure 1: Action density of SU(2) nontrivial-holonomy calorons with increasing 
'instanton radius' (left to right) at fixed temperature. Figures are taken from a 
paper by Kraan and van BaaL The peaks of the action density coincide with the 
positions of constituent BPS monopoles. 

determined by the value of y44(|x| — >■ oo. In the following we only consider SU(2) 
nontrivial-holonomy calorons with no net magnetic charge. Analytical expression 
for SU(2) caloron solutions of trivial (nontrivial) holonomy can be found in [23 
(I2H1 Oil 123), see also Fig. 1.. Since calorons are BPS saturated or self-dual their 
energy-momentum tensor vanishes identically. 

An SU(2) caloron of topological charge k has a classical Euclidean action Scai,2 = 
For trivial holonomy the one-loop effective action of a charge-one caloron is 
given as [HO] 



where p and T denote the 'instanton' radius and temperature, respectively. For 
large g and small enough p the trivial-holonomy caloron thus sizably contributes to 
the partition function of the theory. The one-loop effective action of a nontrivial- 




(5) 
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holonomy caloron is 

^eff cx TV (6) 

where V denotes the spatial volume of the system. In the thermodynamic limit V — > 
oo nontrivial-holonomy calorons thus do not contribute to the partition function. 
As we will show below, however, the thermodynamic limit is not physical due to 
a domanization of the ground state of the theory. The suppression of nontrivial- 
holonomy calorons in the partition function is thus governed by the size of a typical 
domain. 

It was shown in that SU(2) calorons are unstable under one- loop quantum 
fluctuations. Namely, for a holonomy close to trivial there is an attractive potential 
between constituent monopole and antimonopole. In the opposite case the potential 
is repulsive implying the dissociation of the caloron into a monopole-antimonopole 
pair. In a mesoscopic level, an isolated monopole arises at a point in space where 
four or more Higgs-field domains meet [H^ . 

Since the action density of a caloron is T-dependent the action density of the 
macroscopic, adjoint Higgs field (p should be T-dependent through the T-dependence 
of the configuration (j). 

The effective theory describing the (electric) phase macroscopically is an adjoint 
Higgs model: 

In Eq. ((7j) Ve denotes the potential responsible for the stabilization of 0. The covari- 

13 



ant derivative is defined as Vpcj) = dp + ie[(f), Op], the field strength as G^^ = G^j^t"-, 
where = (9^a" — d^al^^ — ef"'''^a''^al, e denotes the effective gauge couphng con- 
stant, and trNt^t'' = 1/2 5"^. While the effect of nontrivial topology is described 
by the scalar sector of the effective theory the curvature G^^, is generated by the 
topologically trivial fluctuations. 

For future work we propose the following nonlocal definition for the phase 
of a given SU(2) block (p. 



ti2 I dp d'x A" F^,((r, 0)) [(r, 0), (r, x)]F^,((r, x))[(r, x), (r, 0)] + ■ ■ ■ . (8) 

The dots in Eq. (jHJ denote the contributions of higher n-point functions, and the ~ 
sign indicates that this expansion very likely is asymptotic at best as a powers series 
in a dimensionless parameter ^. This, however, is not an obstacle to determining 
(f)'s phase and modulus [33]. Each block receives a nontrivial phase by the cor- 
responding SU(2)-embedded trivial-holonomy caloron (or anticaloron A^) over 
which the correlator in Eq. (jH)) is evaluated^. In the definition Eq. (jH} [(r, 0), (r, x)] 
denotes a Wilson line in the fundamental representation which is taken to be along 
a straight path connecting the two points (t, 0) and (r, x): 



[(r,0),(r,f)]=Pexp 



/.(r,x) 

i / dypAp{y) 

J{t,0) 



(9) 



In a lattice simulation at finite temperature T the average in Eq. (jHJ can be computed 
by using an ensemble of cooled configurations whose action is an integer multiple of 



The topological charges of or may, in principle, vary from block to block. 
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87r^/(7^^. Local gauge-singlet composites such as the gluon condensate 



(10) 



are thermodynamically irrelevant for the following reasons: Since they do not cou- 
ple to the topologically trivial sector they do not influence the mass spectrum of 
fluctuations 5ap. Moreover, a singlet composite, arising from noninteracting trivial- 
holonomy calorons, would have zero energy density and pressure because of the BPS 
saturation: a situation which cannot be changed by interactions mediated by the 
trivial sector due to the missing gauge charge. The situation is different though if an 
axial anomaly, arising from integrated-over chiral fermions, is operative. In this case 
the composite in Eq. (jl(J|) determines the mass of the axion, and thus it is visible. 

The key question now is whether the potential Ve in Eq. ((7)) is uniquely de- 
termined by our basic assumption. What are the properties of the field that 
can be deduced? In thermal equilibrium (j) must be periodic in Euclidean time r 
(0 < r < 1/T). Since (f) describes the ground state of the thermal system its modulus 
101 must not depend on t,x but should depend on T. Since cf) is built of noninter- 
acting, self-dual configurations (zero energy density and pressure) it must also be 
pressure - and energy- free. This is the case if and only if 0(r) (not its modulus!) is 
BPS saturated, that is, it solves the following equation 



dA = Ve 




^The nontrivial-holonomy part is then cooled away. 
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where ve is a 'square root' of the potential Ve'- 

VE{4>)=ir^v^EVE. (12) 

The above properties fix the potential uniquely to be VE(0)trN = A|;trN(0^)~^ As it 
turns out, a (winding) solution to Eq. (fTT|) is quantum mechanically and thermody- 
namically inert and thus can be used as a background to the macroscopic equation 
of motion for the trivial-topology sector of the theory. 
The equation of motion 

V^G^, = 2ie[<P,V,<P\, (13) 

which follows from the effective action ((Tj), determines a configuration a^^. For 
a^^ to describe the ground state of the theory it needs to be pure gauge, that is, 
G^,p[a!'^] = 0. Otherwise the invariance of the thermal system under spatial rotations 
would be spontaneously broken. It will turn out that such a pure-gauge solution 
— exists for = 0. As a consequence, the action density in Eq. ((7)) when 
evaluated on 0, a^^ reduces to the potential Ve- We thus describe on a macroscopic 
level interactions between trivial-holonomy calorons as mediated by the topologi- 
cally trivial sector. Namely, the vanishing ground-state energy density (pressure) 
of noninteracting trivial-holonomy calorons is shifted from zero to Ve{4>) (— Ve(0)). 
Moreover, a macroscopic holonomy arises which indicates that (unstable) nontrivial- 
holonomy calorons are generated by gluon exchange and, as a consequence, that 
isolated magnetic monopoles occur. This precludes our conceptual discussion of the 
ground-state physics. 
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An adjoint Higgs field breaks the SU(N) gauge symmetry down to U(l)'^'^ at 
most. Whether SU(N) gauge symmetry is broken maximally or submaximally is 
decided by a boundary condition to the BPS equation (jlip set at an asymptotically 
high temperature. Interacting calorons emit and absorb gauge-field fiuctuations 
5 a p. To discuss their quasiparticle mass spectrum a gauge transformation to a^^ = 
(unitary gauge) must be performed. We will show explicitly for the SU(2) case 
that such a transformation is on the one hand nonperiodic but on the other hand a 
symmetry transformation for all thermodynamical quantities. This is true since the 
transformation does not affect the periodicity of the fluctuations 5ap (no Hosotani 
mechanism [Sn])- The nonperiodic gauge transformation maps the Polyakov loop 
from —1 to 1, therefore generates a global electric center transformation and thus 
interpolates between the two physically equivalent ground states of the theory. As 
a consequence, the global symmetry Zj^oicc is spontaneously broken and hence the 
electric phase is deconflning. The generalization to arbitrary N is straight forward. 

There are tree-level heavy (TLH) and tree-level massless (TLM) modes in 5ap. 
Due to the T dependence of on-shell TLH modes are thermal quasiparticles. 

Due to the T dependent Higgs mechanism and the T dependent ground-state 
energy, which are both generated by the macroscopic fleld 0, implicit temperature 
dependences arise in a loop expansion of thermodynamical quantities. 

To guarantee in the effective theory the validity of the Legendre transformations 
between thermodynamical quantities, as they can be derived from the partition 
function of the underlying theory, thermodynamical self-consistency has to be de- 
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manded. This condition determines the temperature evolution of the effective gauge 
couphng constant e with temperature. As we will see, there is an attractor to this 
evolution which is the constancy of e except for a logarithmic pole at a temperature 
Tc^E- We thus recover in the effective theory the ultraviolet-infrared decoupling that 
follows from the renormalizability of the underlying theory. 

The approach to the ground-state dynamics is an inductive one. Namely, we 
first define a potential Ve{(P) and subsequently show that this definition implies the 
above properties of (j), a small action for calorons at the temperature where they are 
assumed to first form the field 0, and the existence of a pure-gauge solution a^^ to 
Eq. (|13|). The thermal system decomposes into a ground state and a part represented 
by very weakly interacting quasiparticle fluctuations ^. One can consider the former 
as a heat bath for the latter at low temperatures and vice versa at high temperatures. 

2.2 Caloron 'condensate', macroscopically 
2.2.1 The case of even N: Ground-state physics 

We first address the macroscopic dynamics of the adjoint Higgs field when N is 

even. The case of odd N is discussed in Sec. 12. 2.^ We may always work in a gauge 
^We compute two- loop correction to the pressure in |37| . 
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where is SU(2) block diagonal: 





02 

■•• 



(14) 



V ■ ■ / 

In Eq. (IT^ each field 0^, (/ = 1, ■ ■ ■ , N/2), lives in an independent SU(2) subalgebra 
of SU(N), and we define the SU(2) modulus as 



7 ,2 1 7 2 
5z| = -tr2 0i , 



(15) 



The potential Ve in Eq. ((7)) is defined as 



Ve = tTt^vl^VE = A| tlN 



(16) 



where A^; is a fixed mass scale generated by dimensional transmutation. It is im- 
portant to note already at this point that there is only one independent mass scale 
describing the thermodynamics in all phases of the theory. 
We define f^; as follows: 



Ve = ^A| 



Al0l/|01 











Ai02/|02|' 





(17) 



V ■ ■ / 

where Aj, {i = 1,2,3), denote the Pauli matrices. This definition is modulo global 
SU(2)-block gauge transformations. This global symmetry (the 'direction' of wind- 
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ing along a U(l) circle around the group manifold od SU(2)) will translate into 
a gauge symmetry once the theory condenses magnetic monopoles, see Sec. l3.2l 
In SU(2) decomposition the solution to the BPS equation (fTT|) reads 



A3 

^ \3exp{-27iiTK{l)XiT) {U 



\27rT\K{l)\ 

where K{1) is a non-zero integer. The solution in Eq. ()18j) is periodic in r and 
depends on T. The set {K{1), . . . , K(N /2)}, which is a boundary condition to 
Eq. |TT|) at the large temperature Tp, determines the value of the potential at a 
given temperature. It also specifies to what extent the SU(N) gauge symmetry is 
spontaneously broken by caloron 'condensation'. For example, the sets {1, 1} and 
{1,2} break SU(4) down to SU(2) xSU(2) xU(l) and V{lf, respectively. Out of 
15 gauge-field modes 7 modes remain massless in the former and 3 modes in the 
latter case. For a description in terms of a given SU(N) Yang-Mills theory the set 
{i^'(l), . . . , i^'(N/2)} has to be measured, see also the discussion in Sec. 12. 51 For 
definiteness and simplicity we assume in the following that the gauge symmetry 
breaking is maximal in such a way that the potential Ve{4>) is minimal (MGSB). 
This corresponds to the boundary condition {K{1), . . . , K(N/2)} = {1,2,..., N/2} 
or a (local) permutation thereof. 

Let us now verify that the solution in Eq. p8|) is quantum mechanically and 
thermodynamically stabilized. Assuming MGSB, the following ratios are obtained 

T2 |0;|2 ^ ^ ^ 

where the dimensionless temperature is defined as A^; = 2'kT / h^E- For N not too 
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large we have A^; ^ 1, see Sec . 12 .41 As one can infer from Eq. ()19|) . the mass mf = 
f^l^^^l^B of collective caloron fluctuations is much larger than T and the compositeness 
scale \(f)i\. The off-shellness of quantum fluctuations of the field (pi is cut off at this 
scale in Minkowskian or Euclidean signature as 

- mf\ < \4>iW or pl + m^<\4)i\\ {pI > 0) . (20) 

So if no off-shellness in Minkowskian or Euclidean signature is allowed for on the 
one hand. On the other hand, statistical fluctuations of on-shell 0-particles are 
strongly Boltzmann suppressed and thus negligible. We conclude that the solution 
(j) in Eq. (fTSj) is stabilized against fluctuations 6(f) and the potential Ve in Eq. p6|) is 
a truly effective one. Thus (p is nothing but a background for the thermodynamics of 
the topologically trivial sector of the theory. As we will see in Sec. 12. 41 topologically 
trivial quantum fluctuations 6ap generate only a tiny correction to the tree-level 
value Ve{(I)). 

Before we investigate the properties of the fluctuations 5ap let us complete our 
construction of the ground state. The ground-state configurations a^^ needs to be 
a pure-gauge solution to the classical equation of motion Eq. (jl3j) . In order for a^^ 
not to break the rotational invariance of the system it needs to be pure gauge. 
Inserting the background (|18|) (winding numbers: {/ = 1, . . . , / = N/2} for MGSB) 
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into Eq. (fT^. we obtain the following pure-gauge solution: 

/ \ 
Ai •■■ 



2Ai 




(21) 



Moreover, we have 



Pp0 = O 



(22) 



on 0, a^^. A remarkable thing has happened: On a macroscopic level we describe 
the generation of a nontrivial holonomy by interactions between trivial holonomy 
calorons, mediated by trivial-topology fluctuations, in terms of a macroscopic holon- 
omy associated with a^^! For the microscopic physics this implies the generation of 
(rare) nontrivial-holonomy calorons and their subsequent dissociation into magnetic 
monopoles. On a mesoscopic level, this is nothing but the Kibble mechanism for 



monopole creation (HB] arising from the domanization of color orientations of 0. 
Moreover, the vanishing pressure and energy density of a hypothetical ground state 
composed of noninteracting trivial- holonomy calorons, is shifted to =F^£;(0) with 



V^i^(0) = f A|rN(N+2) 



(23) 



by gluon exchange (insert Eqs. (PT|) and into Eq. ((7))). 

Let us now split the topologically trivial part in Eq. further into the ground- 
state part a^^ and fluctuations 6ap: 



(24) 
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To make the mass spectrum of the fluctuations 6ap visible it would be desirable to 
work in unitary gauge where a^^ = and thus no coupling of 6ap to the background 
takes place. 

The gauge rotation = e*^, which transforms (p and ap according to 



ap apn + -{dpn^)n 



(25) 



from winding gauge to unitary gauge is for MGSB given as 



— ttAiTt 



















-2'kXiTt 







01 

02 

■•• 



(26) 



\ ■ ■ / V • • / 

Notice that the gauge transfromation VL as parametrized by Eq. fj26|) is not periodic 
due to its first, third, fifth, ... block being antiperiodic in r. Is a nonperiodic gauge 
transformation physically admissible? Let us discuss this for the SU(2) case only. 
We can make = exp [ ] periodic at the expense of sacrificing its smoothness 



at the point r = by 



where Z{t) is a local (electric) Z2 transformation of the form 



(27) 



Zir) = 2Q{r-—)-l 
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and B denotes the Heavyside step function: 

0, (x<0), 

e(a;) = <! i, (x = 0), ■ (29) 

1 , (x > 0) . 

Applying Q! to = a^^ + 6ap, where 6ap is a periodic fluctuation in winding gauge, 
we have 



ap Q^al^ + 5ap)Q+^-dpQ^Q 



Q\a'f + 6ap)Q + - ({dpQ^)Q + {dpZ{T))Z{T 



= n^apn + ^6{T-^)Z{r) 

= n^dapVl . (30) 

Since ^"''(0) = — fi"l'(^) = fl{0) = — fi(^) we conclude that if the fluctuation 6ap is 
periodic in winding gauge it is also periodic in unitary gauge. It thus is irrelevant 
whether we integrate out the fluctuations 6ap in winding or unitary gauge in a loop 
expansion of thermodynamical quantities: Hosotani's mechanism |36| does not take 
place. What changes though is the Polyakov loop evaluated on the background 
configuration a^^: 

P[aJ^] = -1 ^ P[0] = 1 . (31) 

We conclude that the theory has two equivalent ground states and that the global 
electric ^'a.dec symmetry is spontaneously broken. We thus have shown that the 
elecric phase is, indeed, deconfining. The generalization of this result to SU(N) with 
N even is straight forward. 
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In unitary gauge we have 



Mr) 




A3- 



(32) 



Thus the field is constant and diagonaL Moreover, we have 



(33) 



The gauge- covariant kinetic term for in the action Eq. (|7j) reduces to a sum over 
mass terms for the TLH modes contained in 6ap. The TLH (TLM) modes are 
massive (massless) quasiparticles associated with three (two) polarization states. As 
we will show in by computing the two-loop correction to the thermodynamical 
pressure for N=2 these quasiparticles are practically noninteracting for sufficiently 
large temperatures. 

2.2.2 The case of odd N 

If N is odd then a decomposition of (p into SU(2) blocks only is no longer possible. 
One of the SU(2) blocks in Eq. (jl4p is then replaced by an SU(3) block. Impos- 
ing MGSB and counting the number of independent stable and unstable magnetic 
monopoles microscopically on the one hand and macroscopically in the electric and 
the magnetic phase on the other hand, see Sec. (j2.6j) . we conclude that temporal 
winding should take place within the SU(3) block as in Eq. (fTHj) but now in each 
of the two independent SU(2) subalgebras only for half the time. The first SU(2) 
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subgroup is generated by Aj, (i = 1, ■ ■ ■ , 3) where 



Ai 







1 







A, 



-i 



i 



A, 







1 



0-1 



(34) 



The second SU(2) subgroup is generated by A,, {i = 1, - ■ ■ , 3) where 



Ai 



1 







yl Oj 



I 



A, 



-i 




y 2 y 



A, 



1 







(35) 



The solution of the BPS equation (fTT|) for the SU(3) block can be obtained by 
applying the following prescription to a single block (fT8|l of odd winding number 
Kipj in Eq. p4j] : generate two configurations and (^1^2 by replacing the Pauli 
matrix Aj in Eq. pSj) by Aj in the first half period, < r < 1/(2T), and by zero in 
the second half period, 1/(2T) < r < 1/T, and by replacing the Pauli matrix Aj in 
Eq. ffTHj) by zero in the first half period, < r < 1/(2T), and by Aj in the second 
half period, 1/(2T) < r < 1/T. Add and ^1^2 to generate a new solution to 
the BPS equation (fTT|) over the entire period. Within this block a member of the 
third, dependent SU(2) algebra is generated at r = 0, 1/(2T). To which block I of 
odd winding number K{t) this prescription is applied is a boundary condition to the 
BPS equation. For simplicity we will proceed in this paper by only quoting results 
for the case N=3. 

It is clear that for even > 6 a decomposition of may contain an even 
number of SU(3) blocks besides SU(2) blocks. For definiteness, we only consider the 
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decomposition into SU(2) blocks as proposed in Eq. (fT^ . 



2.3 Tree-level mass spectrum of TLH modes 

In absence of radiative corrections only N(N— 1) TLH modes acquire mass by the 
adjoint Higgs mechanism. TLM modes aquire tiny screening masses radiatively. 
This effect will be discussed in Sec. l2.7l 

In unitary gauge, the TLH mass spectrum calculates as 



mi = -2e' tr [0, t'] [0, , (k = !,■■■, N(N-l)) . 



(36) 



The SU(N) generators t*^, which are associated with the TLH modes, are 



(/ = 1,---,N; J>/; r,s = l,---,N) 



(37) 



By virtue of Eqs. ^ we obtain: 



Al 



mjj = mij = e (0/-0j) = e ^ { 



1/2 J J 12 ) 



//2 V(j+i)/2; ' 



(/even, J even) 

(/odd, J odd) 

(/odd, J even) 

(/even, J odd) . 

(38) 



For N=3 we have 



12 



m;^3 = m 



2 

12 



A3 

27rT 
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= 4e 



A3 

2 

27rT 



or 
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^12 = ^23 = ^12 = "^23 = e 



2 AL 
27rT 



mis = ^i3=4e'^. (39) 



2.4 Thermodynamical self-consistency and e{T) at one loop 

The TLH modes (5a^ are thermal quasiparticle fluctuations on tree-level since their 
masses mjj and rrijj, given in Eqs. (|HH|) and (jHUj) . are T dependent. Moreover, 
the ground-state pressure, given by — Ve, is linearly dependent on T, see Eq. (j2SI)- 
Thermodynamical quantities such as the pressure, the energy density, or the entropy 
density are interrelated by Legendre transformations as derived from the partition 
function associated with the underlying SU(N) Yang-Mills Lagrangian. To assure 
that the same Legendre transformations are valid in the effective electric theory, 
where ground-state pressure and particle masses are temperature dependent, a con- 
dition for thermodynamic self-consistency needs to be imposed. In general, this 
condition assures that T-derivatives of quantities that enter the action density of 
the effective theory (in our case the TLH masses and the ground-state pressure) 
cancel one another. 

Let us formulate this condition at one-loop accuracy. It is convenient to work 
with dimensionless quantities. The quantity (mass over temperature) is defined 
as 



ak = Cka where « = , (/c = 1, ■ ■ ■ , N(N - 1)) , (40) 

and the coefficient is equal to the square root of one of the numbers appearing 
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to the right of the curly bracket in Eq. (j5H|) or, for N=3, to 1 or 2, compare with 
Eq. RecaUing our definition of a dimensionless temperature 



Ak = 



2'kT 



(41) 



we have 



a = 27reA 



-3/2 



(42) 



Let us also define the (negative definite) function P{a) as 



_ roc 

P(a) = dxx^ log[l — exp(— V + a?)] . 
Jo 



(43) 



Ignoring higher loop corrections, the total thermodynamical pressure P{Xe) associ- 
ated with the diagrams in Fig. Q calculates as 



N{N-1) 

2(N-l)P(0) + 3 

k=l 



Ar /N 



(44) 



In Eq. (j44p we have neglected the 'nonthermal' contribution —AVe to the one-loop 
bubbles in Fig. El which is estimated as (quantum fluctuations are cut off at the 
compositeness scale |0|) 

Q ^1^1 / ^ \ 1^14 

(45) 



1287r2 • 



To neglect —AVe is well justified for sufficiently small N, say N<20, since we have 
AVe 



Ve 



< - 1)^ f MV ^ (N^ - 1)^X-/ , 



(46) 



and we will show later that the minimal temperature in the electric phase X^^e is 
much larger than unity. 



29 




Figure 2: Diagrams contributing to the pressure when radiative corrections are 
ignored. A thick hne corresponds to TLH modes and a thin one to TLM modes. 
The cross depicts the ground-state contribution arising from caloron 'condensation'. 



For N=3 we obtain 



P{Xe) = -A^I^ [4P(0) + 3 (4P(a) + 2P(2a))] + 2A^ | (N=3) . (47) 

A particular Legendre transformation following from the partition function of the 
underlying theory maps pressure into energy density as 

dP , , 

P-T--P. (48) 

For Eq. ()48|1 to hold also in the effective electric theory the following situation has 
to be arranged for. only the explicit T-dependence in P, arising from the explicit 
T-dependence of the Boltzmann weight, should contribute to the derivative ^ while 
implicit T dependences of gauge-boson masses and the ground-state pressure ought 
to cancel one other. This condition is expressed as 



daP = 0. (49) 

Before we proceed let us recall that the A^; dependence of the ground-state pressure, 
as indicated in Eqs. (j^^ and (HH), could be expressed in terms of a dependence 
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Figure 3: The function aD{a). 



on the mass parameter a by virtue of Eq. (jlzj) if the A^; dependence of the gauge 
couphng constant e was known. Keeping this in mind, we derive from Eqs. 
(jH}, and ipTjl the following evolution equation 



N{N-1) 



For N=3 we have 



daXi 



^^^^''(Z}(a) + 2Z}(2a)) , (N=3) 



(50) 



(27r) 



6 



(51) 



The function D(a) is defined as 



D(a) = 



Jo 



dx 



X 



\/x'^ + a'^ exp{y/x'^ + a^) — 1 



(52) 



Eqs. (Hljl and (|KT|) describe the evolution of temperature as a function of tree-level 
gauge boson mass. The right-hand sides of these equations are negative definite since 
the function D{a) in Eq. ()52j) is positive definite, see Fig.El As a consequence, the 
solutions Xe{cl) to Eqs. (|3U|) and (jHH) can be inverted to ^(A^;). In Fig. HI a solution 
for N=2 subject to the initial condition Xe,p = A£;(a = 0) = 10'^ is shown. We 

31 



have noticed numerically that the low-temperature behavior of Xe{0') is practically 
independent of the value Xe,p as long as Xe,p is sufficiently large. Let us show this 
analytically. For a sufficiently smaller than unity we may expand the right-hand 
side of Eq. (j5(J|) only taking the linear term in a into account. The inverse of the 
solution is then of the following form 



The dependence in Eq. thus is an attractor. So at whatever asymptotically 
high temperature the formation of the adjoint Higgs field (j) out of noninteracting 
trivial-holonomy calorons is assumed does not influence the behavior of the theory 
at much lower temperatures. This result is reminiscent of the ultraviolet-infrared 
decoupling property of the renormalizable, underlying theory. 

Notice that Fig. Q and Eq. (I^Uj) imply that there are fixed points of the evolution 
\e{(i) at a = and a = oo. The points Xe^p = Xe^ci = 0) and Xe,c = Xe^o. = oo) 
are associated with the highest and the lowest attainable temperaturesin the electric 
phase, respectively. In a bottom-up evolution no information can be obtained about 
Xe,p if the temperature that is maximally reached is sufficiently smaller than Xe,p- 
In a top-down evolution, where Xe,p is set as a boundary value, the prediction of 
Xe^c is independent of Xe^p- These two statements are immediate consequences of 
the existence of an attractor in the thermodynamical evolution. 




(53) 



For Xe sufficiently smaller than Xe^p this can be approximated as 



a oc A^;' 



(54) 
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Numerically inverting the solution Xe{cl) to ^(A^;), the evolution of the gauge 
coupling constant e can be computed using Eq. 



We show the result in Fig.Elfor N=2,3. Before we interprete this result a remark 
on the interpretation of the effective gauge coupling constant e for T ~ Tp is in 
order. Since e determines the strength of the interaction between nontopological 
gauge field fluctuations Sap and the coherent caloron state there is, in general, no 
reason for it to be equal to the gauge coupling constant g of the fundamental Yang- 
Mills theory. However, for temperatures very close to Tp, where (j) is assumed to 
form (see also Sec. l2.5p . the coupling constant e should be roughly equal to g. From 
Fig. El we see that the effective gauge coupling constant e ~ ^ evolves to values larger 
than unity shortly below the initial temperature Xe,p- This is in agreement with 
our assumption that trivial-holonmy SU(2) calorons (of sufficiently small 'instanton' 
radius) have a large action and thus contribute sizably to the partion function of 
the underlying theory, compare with Eq. (0). For a grandly unifying SU(N) gauge 
theory we argue in Sec. 12. 51 that Xe,p ~ ^ where Mp ^ Ae denotes the cutoff 
scale for the description in terms of a local, four-dimensional field theory. 

There is no handle on the relation between the fundamental gauge coupling g 
and e after the condensate has formed and is sustained by interactions between the 
trivial-holonomy calorons. Since collective effects are strong they can make up for 
a large value of the action of an isolated caloron generated by a small value of g. 




(55) 



33 



^ 200 
150 
100 
50 



0.5 1 1.5 2 2.5 3 3.5 

a 

Figure 4: The solution A£;(a) to Eq. (48) for N=2 subject to the boundary condition 
Xe^p = \E{a = 0) = lO^. 

This may justify the use of universal perturbative expressions for the beta function 
in lattice simulation at low temperatures. 

Notice that the dependence of a on A^; in Eq. (jMj) is canceled in the depen- 
dence of e on As such that a plateau is reached quickly in Eq. We interprete 
the fact that the gauge coupling constant e remains constant for a large range of 
temperatures as another indication for the existence of spatially isolated and con- 
served magnetic charges in the system, see also Sec. 12.61 During the relaxation of 
e to its plateau value constituent BPS magnetic monopoles residing in dissociat- 
ing nontrivial-holonomy calorons form as isolated defects [211 I2H1 1^2] • Since the 
interaction between a monopole and an antimonopole, as mediated by TLM modes, 
is screened these defects need not be considered explicitly in the effective, 

thermal theory discussed in the present work. Implicitly, their presence is accounted 
for here by the holonomy of the background field af'^ . 

The effective gauge coupling constant e runs into a logarithmic (needle) pole at 
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Figure 5: The low-temperature evolution of the gauge coupling e in the electric 
phase for N=2 (grey line) and N=3 (black line). The gauge coupling diverges log- 
arithmically, e oc — \og{XE — K,e), at \e,c = 11-65 (N=2) and \e,c = 8.08 (N=3). 
The plateau values are e = 5.1 (N=2) and e = 4.2 (N=3). 

Xe,c of the form 

c^Xe) oc - log(Aij - Xe,c) , (56) 

compare with Fig.^J 

The plateau values for e are e ~ 5.1 and e ~ 4.2 for N=2 and N=3, respectively. 
For a given N they do not depend on where the boundary condition Xe,p = Xe{cl = 0) 
is set if Xe,p is sufficiently larger than Xe,c- For N=2 we have Xe,c = 11-65 and for 
N=3 we have Xe,c = 8.08. It thus is self-consistently justified to neglect the one-loop 
quantum corrections to V£;(0) as they are estimated in Eq. Naively, one would 
conclude that the large plateau values would render the fiuctuations 6ap to be very 
strongly coupled and that radiative corrections to the thermodynamical potentials 
would thus be uncontrolled. This, however, does not happen due to the fact that 
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the TLH modes acquire masses by the Higgs mechanism which are proportional to 
e and due to the existence of a compsiteness scale |0|. The latter constrains the 
momentum p of quantum fluctuations in Sap as 

\p'^-ml\ < (TLH, Minkowskian) , < (TLM, Minkowskian) , 

pl + ml < (TLH, Euclidean) , < (TLM, Euclidean) . (57) 

Since in nonlocal two-loop contributions, see Fig.[3 each TLH line is on-shell for 
e > 1 the effect of a strongly coupled vertex is compensated by the very small phase 
space that is allowed for the progagation of a TLM mode coupling to the TLH 
mode. The center-of-mass energy flowing into or out of a four-vertex is constrained 
in addition to Eq. ()57|) to be smaller than We show in that the two-loop 
contributions to the pressure for SU(2) are, depending on temperature, at most 
~ 0.1% of the one- loop result. 

2.5 What is Tp? 

At temperatures larger than the highest attainable temperature Tp in the electric 
phase (corresponding to a = 0) a grandly unifying SU(N) gauge symmetry, which 
generates all matter and its (nongravitational) interactions in the Universe at lower 
temperatures, would be unbroken. We assume here that gravity is a perfectly clas- 
sical theory up to the Planck mass Mp ^. The perturbative phase of the SU(N) 

Yang-Mills theory at T > Tp would have a trivial vacuum state represented by 
^This assumption is usually made in field-theory models of cosmological inflation. 
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weakly interacting quantum fluctuations. The momenta associated with these fluc- 
tuations can be maximally as hard as some cutoff scale at which the four- dimensional 
setup ceases to be reliable. Common belief is that this cutoff scale is Mp. 

The highest temperature T^^^^^ attainable in the perturbative phase thus is com- 
parable to Mp, T^utoff ~ ^p- At temperatures rnaging between T^utoff and Tp per- 
turbative vacuum fluctuations would generate a cosmological constant A^^.^^ given 
as 

A.„._ ~ M| . (58) 

At T^utoff the vacuum energy density Mp would be comparable to the thermal energy 
density of on-shell fluctuations ~ ^cttofi- While the former is constant the later dies 
off quickly as the Universe cools down. Thus for T slightly smaller than T^utoff the 
vacuum energy density dominate the expansion hence the Universe would rapidly 
decrease its temperature as 

T~exp[-MpAt]T,,,„ff (59) 

where At = t—tp and tp ~ Mp^. A sudden termination of this Planck-scale inflation 
would occur at T = Tp where the field comes into existence. While the radiation 
component of the total energy density is continuous across the phase boundary at 
Tp the energy density of the ground-state would be discontinuously reduced from 
Mj, to ~ TpAl-j^^, see Eq.^. On the one hand, this release of latent heat is a 
characteristic for a (strong if Tp <C T^utos) 1"* order transition. On the other hand, 
the order parameter a for the onset of the electric phase is continuous, see Fig. © 
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(screening masses are comparable on both sides of the phase boundary, see Eq. (f7T|) ). 
But this is signaUing a 2""^ order phase transition. There is only one way to avoid 
this contradiction: The phase boundary at Tp needs to be hidden beyond the point 
T,,t„ff ~ Mp, that is, Tp > T,^,,g. 

The reader may object that our conclusion about the 2"'* order of the phase 
transition (order parameter a) is resting on a one-loop analysis of the gauge- coupling 
evolution. Usually, it is understood that such a mean-field treatment breaks down 
close to a 2"^* order transition due to fact that long range correlations mediated 
by low- momentum quantum fluctuations become important. In the electric phase 
these long-range correlations are, however, contained in the field which does not 
fluctuate at any temperature Te,c < T < Tp. Due to |0| being a cutoff for the 
quantum fluctuations of the TLM and TLH modes and due to the fact that |0| dies 
off as T- the Ion g-range correlating effects of TLM and TLH quantum fluctuations 
can safely be neglected if Tp ^ Ae- The above discussion and conclusion thus are 
valid. 

So far we had in mind the simplified case of MGSB at Tp in grandly unify- 
ing SU(N) Yang-Mills theory. We know from experiment, however, that gauge- 
symmetry breakdown at Tp is not maximal in Nature. If the gauge-symmetry 
breaking by (p ^'^ Tp is submaximal, however, the same contradiction between the 
orders of the phase transition arises for Tp < Mp. In this case a fundamental gauge 
symmetry SU(N) would be broken to a product of group with factors SU(M) M<N 
or U(l). Recall, that submaximal gauge-symmetry breaking in the electric phase 
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takes place if SU(2) blocks in with equal winding number are generated at Tp, see 
Eq. (dH). 

The theory would then condense magnetic SU(M) color and magnetic U(l) 
monopoles at the temperature Te^c- While condensates of the latter are described by 
complex scalar fields, see Sec. l3.H condensates of the former are, again, described by 
adjoint Higgs fields. Maximal or submaximal breakings of the residual SU(M) gauge 
symmetries would be possible at the electric-magnetic transition of the fundamental 
SU(N) theory. For the effective description of SU(M) thermodynamics at T < Te^c 
this boundary condition effectively is set during the phase transition at Te^c and 
not at Tp. By matching the thermodynamical pressures at Te^c the scale A'^. of the 
effective theory SU(M) is determined in terms of the scale A^; of the fundamental 
SU(N) theory and the pattern of symmetry breaking at Tp and Te,c- After a se- 
quence of such matching procedures has taken place (in which residual U(l) factors 
have thermodynamically decoupled) a hierarchy between A^; and the scale A'e' of 
an effective SU(L) theory {L <^ N), seen experimentally at low energies (or local 
temperatures for that matter), is generated. 

To summarize, we provided an argument that in a grandly unifying and four- 
dimensional SU(N) Yang- Mills theory the dynamical generation of the adjoint back- 
ground field must take place at a temperature T^utoff ~ Mp where the local field- 
theory description breaks down. Any lower SU(L) gauge symmetry, which is gen- 
erated from the SU(N) theory by a sequence of condensations of (color) monopoles 
and confining transitions, is matched to its 'predessesor' theory in 2'"* order like 
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phase transitions. The scale of this SU(L) gauge theory can be much lower than the 
scale A^; of the fundamental SU(N) theory. 

2.6 Stable und unstable magnetic monopoles 

Due to the presence of an adjoint Higgs field in the electric phase there are 't Hooft- 
Polyakov magnetic monopoles jSllH] which are centered at the isolated zeros of 
On a mesoscopic level, these zeros occur at points in space where four or more color- 
orientation domains of a given block 0; meet [23] • Microscopically, BPS monopoles 
|in] are contained within decaying nontrivial-holonomy calorons ^. 

Close to Tp we have ^ ~ e, and the following processes take place: Trivial- 
holonomy calorons grow rapidly in size, start to overlap, and thus generate calorons 
with holonomy by their interactions. If this holonomy is sufficiently large then 
two following processes take place: (i) SU(2) nontrivial-holonomy calorons of the 
same embedding in SU(N) decay independently into their constituent magnetic 
monopoles and antimonopoles, and (ii) SU(2) nontrivial-holonomy calorons gen- 
erated from trivial-holonomy calorons of different SU(2) embeddings^ in SU(N) do 
not decay into constituent monopoles and antimonopoles since they would have to 

live in instable superpositions of the embeddings of the asymptotic trivial-holonomy 

® A perturbative analysis of 1-loop radiative corrections to an isolated instanton were performed 

in |47| . In |82j this was done for the nontrivial-holonomy caloron. As a result a repulsive potential 

for the constituent monopole and antimonopole was obtained for a sufficiently large holonomy. 
^Recall, that we have assumed that these calorons come with different topological charge to 

obtain maximal symmetry breaking. 
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calorons. While the former process generates stable magnetic dipoles the latter 
generates instable monopoles and antimonopoles. 

In our macroscopic approach it is hard to see how the size of a typical trivial- 
holonomy caloron changes with temperature after e has reached its plateau value 
since e plays a different role than the fundamental coupling constant g. We may, 
however, infer from lattice simulations that trivial-holonomy calorons are large 
enough to not generate a topological susceptibility on lattices of presently feasible 
sizes. 

A magnetic monopole-antimonopole pair, which is connected by a magnetic flux 
line, becomes a stable dipole if the pair has a sufficiently large spatial separation. In 
this case the monopole and the antimonopole are practically noninteracting jl^I and 
thus are stable defects. If a single monopole is produced then it is unstable unless 
it connects with its antimonopole, produced in an independent collision. N— 1 inde- 
pendent SU(2) subgroups exist and so N— 1 independent magnetic monopoles may 
occur in the case of MGSB. Since the monopole constituents in a caloron are BPS 
saturated we also expect an isolated monopole in a stable monopole-antimonopole 
pair to be BPS saturated. The analytical expression for an SU(2) charge-one BPS 
monopole in a gauge where the Higgs field (f) winds around the group manifold at 
spatial infinity [46j is given as 

1 - K(r) H(r) 
er er 

where r = y/W, and f is a spatial unit vector. The form of the functions K{r) and 
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H{r) is 

Cr 

K{r) = — — — , H{r) = Cr coth(Cr) - 1 . (61) 

sinn(6r) 

In Eqs. the mass scale C is proportional to the asymptotic Higgs modulus 
\4>{\x\ oo)| and the gauge coupling constant e. The mass of a BPS monopole is 
given as [IHl 

M = ^J^{\x\ ^ oo)\ . (62) 
A dual, abelian field strength G^i, can be defined as jH] 

G,. = ^ - ^,U^,<^)W)c ■ (63) 

The expression in Eq. (j63|) reduces to Gfj^^ = df^al—d^afi^ in unitary gauge cp"" = 
Eq. ()6H|1 defines the field strength of a dual photon which couples to the magnetic 
charge 47r/e of the monopole. Both the gauge dynamics involving only dual photons 
and magnetic monopoles and the entire gauge dynamics in the electric phase are 
blind with respect to the magnetic (local in space) center symmetry Z-f^^^^g (the field 
strength G^^, the field (p and the covariant derivative are invariant under center 
transformations and, as a consequence of Eq. so is the dual field strength Gf^^)- 
However, a local-in-time transformation G 2'n,c1oc may transform the Dirac string 
between a static monopole and a static antimonopole in a given SU(2) embedding 
into a Dirac string belonging to a dipole in a different SU(2) embedding. This does 
not violate the conservation of total magnetic charge and certainly has no effect on 
any gauge invariant quantity. In an effective theory, where monopoles are condensed 
degrees of freedom, a local ^n.cIoc transformation should thus be represented by a 
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local permutation of the fields describing the monopole condensates and the gauge- 
field fluctuations which couple to them. In such an effective theory the action thus 
ought to be invariant under these local permutations. 

How can we see the occurrence of stable and unstable magnetic monopoles in the 
macroscopic, effective theory for the electric phase? In winding gauge the temporal 
winding of the SU(2) block in is complemented by spatial winding at isolated 
points in 3D space. By a large, r dependent gauge transformation the monopole's 
spatially asymptotic SU(2) Higgs field is rotated to spatial constancy and into the 
direction given by the temporal winding of the unperturbed block (pi. Its Dirac 
string rotates as a function of Euclidean time. Since this monopole is stable, a 
correlated antimonopole must exist. We arrive at a dipole rotating about its center 
of mass at an angular frequency 27r/T. This rotation is an artifact of our choice of 
gauge. Rotating the dipole to unitary gauge by the gauge function 9i in Eq. (pUj) . 
we arrive at a (quasi)static dipole. There are isolated coincidence points (CPs) 
in time where the lower right (upper left) corner of the Z"" {{I + 1)"') SU(2) block 
(now (/ = 1, ■ ■ ■,N/2 — 1)) together with the number zero of its right-hand (left- 
hand) neighbour are proportional to the generator A3. Coincidence also takes place 
between the first and last diagonal entry in (f). At a CP, spatial winding may take 
place at isolated points in 3D space. Moreover, coincidence also takes place between 
nonadjacent SU(2) blocks and the first and last diagonal entry in 0. The associated 
monopoles are, however, not independent. The spatial winding associated with 
the additional SU(2) generators 'flashing out' at the CPs corresponds to unstable 
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magnetic monopoles. 

Summing up all independent monopole species, we have: 

N N 

- + - - 1 = N - 1 . (64) 

In the case N=3 the field (f) winds with winding number one in each of the two 
independent SU(2) subalgebras for half the time, see Sec l2.2.21 The match between 
these subalgebras happens at the CPs tcp = 0, 1/(2T) where an element of the 
third, dependent SU(2) algebra is generated, see Sec. l2.2.2l Due to these CPs we 
have unstable monopoles. 

2.7 Outlook on radiative corrections 
2.7.1 Contributions to the TLM self-energy 

Let us now investigate for N=2 and at one loop the simplest contribution to the 
polarization tensor for the TLM mode. A complete investigation of two-loop con- 
tributions to the pressure is the objective of jnZj. We work in unitary gauge 
= diag(0i,02), flp^ = 0. This condition fixes the gauge up to U(l) rotations 
generated by A3. This remaining gauge freedom can be used to gauge the TLM 
mode to transversality: di6aj^^ = (radiation or Coulomb gauge). No ghost fields 
need to be introduced in unitary-Coulomb gauge. 

After an analytical continuation to Minkowskian signature^° the asymptotic 



-"^"For the purpose of the present work we do not need the matrix formulation of the real-time 
propagators. 
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propagator of a free TLM mode is given as 



DlZfC'^T) = -5^,Pl i^- + 2n5{e)ns{\ko\/T) ) (65) 



where 



pT _ pT _ pT _ n, 
-'00 ~ Oi ~ iO ~ ' 



Pi = (66) 



and nsix) = exp[L]-i denotes the Bose distribution function. The analytically con- 
tinued asymptotic propagator of a free TLH mode D^^^ab^'^{k, T) is that of a massive 
vector boson 



7-,TLH,/j,0/ 7 rTi\ T ( k^ky 



+ 2Ti5{e-mlj)nB{\k,\/T) 



k"^ — m'jj 



(67) 

The vertices for the interactions of TLH and TLM modes are the usual ones. In 
unitary- Coulomb gauge the 4D loop integrals over quantum fluctuations are cut 
off at the compositeness scale |0|(T) of the effective theory. Thermal fluctuations, 
associated with 3D loop integrals, are automatically cut off by the distribution 
function ns- Let us now look at the tadpole contribution to Ila=3,fip as shown in 
Fig-El This diagram decomposes into a part for the vacuum fluctuations in the 
loop, which has a T dependence only due to the T dependence of TLH masses, and 
a thermal part. Contracting the Lorentz indices, the former can be calculated as 



,.2 



^ .mif f!(iiiij) 
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V2,b 



p, a=3 



Figure 6: A tadpole contribution to the self-energy of the TLH mode. 



while the thermal part reads 



K=3r = ^,mrL dy 



y'^ + (2e)^ exp 



27fX-^'/\/y^ + (2e)2 



(69) 



It is instructive to perform the weak and strong coupling limits in Eqs. fl68p and 
For e < we obtain 



T-rvaCj/i 



327r2 



1 + 16e^ - 



12 

^0(l--log(2e)) 



(70) 



and for e ^ 1 



(71) 



For e ^ there is no vacuum contribution, IT^lf^^ = 0. The thermal part reads 

(72) 



where Ki{x) denotes a modified Bessel function. The weak coupling result for the 
thermal part in Eq. ()71|) coincides, up to a numerical factor, with the perturbative 
expression for the electric screening (or Debye) mass-squared, as it should. In the 
limit of infinite coupling, which is reached due to the logarithmic pole for Xe \ Xe,c, 
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see Eq. fl56|) . the thermal part in Eq. ()72j) vanishes. This agrees quahtatively with 
results obtained in thermal quasiparticle models fitted to lattice data [221 HOI HI] • 
It was found in these models that the Debye mass vanishes for T \ Te,c A 
large and constant value of e, as it is generated by one-loop evolution (compare with 
Eqs. fl54|) and implies that the approximation leading to Eq. (f?^ breaks down 

for high temperatures. It is, however, clear from Eq. ()69p that the weak coupling 
result at O(e^) in Eq. (f7T|) gives an upper bound on the contribution to the screening 
mass-squared at any temperature and any value of the coupling constant. 

We expect that the situation is similar for the nonlocal one-loop diagrams. As 
for the tadpole correction in the polarization operator of a TLH mode there is 
a contribution oc for strong coupling which arises from the vacuum part with 
the TLM mode in the loop. This contribution is, however, suppressed due to the 
constraint that the center-of-mass energy flowing into or out of the vertex must be 
smaller than |0| |37j . 

2.7.2 Loop expansion of the pressure 

Two-loop diagrams contributing to the pressure in a real-time formulation are indi- 
cated in Fig.[7| We do not compute them here but in for N=2. A general remark 

concerning thermodynamical self-consistency is in order already here. Recall, that 
"'^"'^We foretake at this point that the deconfinement phase transiiton seen on the lattice is the 

electric-magnetic transition at Te^c- We will discuss in Sec. 16. 21 why the lattice is not capable of 

measuring infrared sensitive quantities such as the pressure at temperatures below Te^c- 
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on one-loop level we have obtained an evolution equation from the requirement of 
thermal self-consistency daP = 0. This gave a functional relation between temper- 
ature and mass which could be inverted for all temperatures in the electric phase. 
After the relation Eq. (j55|) between coupling constant e and mass a was exploited 
we obtained a functional dependence of the effective gauge coupling constant e on 
temperature. Equivalently, we could have demanded dgP = since e is the only 
variable parameter of our effective theory for the electric phase. This would have 
directly generated an evolution equation for temperature as a function of e. 

Radiative corrections AP to the pressure have a separate dependence on a and 

e, 



thermodynamical self-consistency by demanding daP = one has to express the 
explicitly appearing e in Eq. ()73p in terms of a by means of Eq. ()55p and distinguish 
temperature dependences arising from a simple rescaling and those arising from the 
T dependent ground-state physics. For SU(2) we have 



The first factor on the right-hand sides of Eq. (j74j) arises from rescaling, so only the 
second factor needs to be differentiated: 




(73) 



where AP is a dimensionless function of its dimensionless arguments. To implement 




(74) 




(75) 
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4(CD^OO) 



Figure 7: Two- loop diagrams contributing to the pressure. Thick hues denote prop- 
agators of TLH modes, thin hues those of TLM modes. 

After solving daP = for the term da^E we obtain a modified right-hand side of 
the evolution equation Eq. ()50p . The inverted solution to this evolution equation 
describes the dependence of mass on temperature or, after applying Eq. the 
dependence of e on temperature when two-loop diagrams for the pressure are taken 
into account. The Euclidean momenta Pe of off-shell fluctuations of the TLH- modes 
are constrained by the condition 

pI < - 4 ^) ' (A; = 1, ■ ■ ■ , N(N - 1)) , (76) 

and by the requirement that the total momentum squared flowing into or our of a 
four- vertex cannot be larger than s ee Eq. ()20|) and Eq. ()40|1 . Since at one loop 
e is smaller than unity for T ~ Tp only, we expect TLH-mode quantum fluctuations 
to be absent at temperatures lower than Tp also at higher-loop accuracy. 
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3 The magnetic phase 



3.1 The electric-magnetic phase transition 

In Sec. l2.6l we have discussed how stable and unstable BPS monopoles are generated 
as isolated objects in the electric phase. For definiteness we have assumed MGSB 
by the adjoint scalar cf). The mass of the N/2 stable BPS monopole species is given 
as in Eq. when replacing 0-^0;, (/ = 1, ■ ■ ■ , N/2). As a consequence of the 
evolution of the gauge coupling c^Xe) following from Eq. (j^Uj) the mass of a stable 
monopole vanishes ai \e = \e,c due the logarithmic pole of e. Stable monopoles do 
not carry any Euclidean action at this point, and thus they condense. 

TLM modes, which couple to isolated monopoles with strength g = become 
dual gauge bosons. They couple to the monopole condensates with a strength g, 
which may now continuously vary with temperature, starting with g = ai Xe = 
Xe,c- The TLH modes of the electric phase decouple kinematically at A^; = Xe,c 
since their masses, oc e|0|, diverge. At the onset of the magnetic phase, where N/2 
species of stable monopoles are condensed, we are thus left with an effective Abelian 
theory of N — 1 dual gauge fields a^,^, {k = 1, ■ ■ ■ ,N — 1) and N/2 condensates 
of stable monopoles described by complex scalar fields ipi , (/ = 1, ■ ■ ■ , N/2). The 
temporal winding of these fields is the same as that of the associated SU(2) blocks 
(pi in the electric phase. For N=3 there are two independent condensates of stable 
monopoles. 

What happens to the N/2— 1 independent unstable monopoles (N even, N > 3? 
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Unstable monopoles are generated by gluon exchanges between trivial-holonomy 
calorons in different SU(2) embeddings. We conclude, that at A^; = Xe,c the intact 
continuous gauge symmetry U(l)^"^ of the electric phase becomes a gauge symmetry 
U(l)^'^ which is spontaneously broken as 

U(1)S' ^ U(l)^'"-^ (N > 3) (77) 

for A^; < Xe,c in the magnetic phase. For = 2,3 the spontaneous breakdown 
of continuous gauge symmetry is maximal. The monopole condensate (pk, which is 
associated with the dual gauge- field fluctuation 5a^^,, is defined as 

y;,, (A; = l,---,N/2), 
^k = \ . (78) 

0, (fc = N/2 + l,---,N-l), 

The local ^'n.cIcc symmetry acts on 5a^^^i^k as a local-in-time permutation (see the 
discussion in Sec. l2.6p 

'^k) (^«^(A:+j(r))mod (N-l)5 <^(fc+i(r))mod (N-l)) , (j e Z) . (79) 

In Eq. fl7^ the integer- valued functions j are piecewise constant on extended regions 
of Euclidean spacetime. The symmetry defined in Eq. (f7^ leaves the ground state of 
the system invariant, and thus the discrete, local symmetry ^N^eiec is unbroken in the 
magnetic phase. As a consequence the global Z^^^iec associated with the Polyakov 
loop as an order parameter is also unbroken, see also Sec. 13. 41 
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3.2 Monopole condensates, macroscopically 

The effective theory describing the magnetic phase is constructed in close analogy 
to the effective theory describing the electric phase. Recall, that we assume MGSB 
in the electric phase. Since the condensation of monopoles is driven by their mass- 
lessness the complex scalar fields ipi, which describe the monopole condensates, are 
energy- and pressure-free in the absence of monopole interactions mediated by dual 
gauge-field fluctuations in the topologically trivial sector of the theory. For the N=2 
case the exponent of the phase of the local field ip is defined as 



In Eq. (j8Uj) the dual field strength G^y is the 't Hooft tensor of Eq. (j63p . the (surface-) 
integal is over a spatial 2-sphere of infinite radius, and the average is over the zero- 
mode deformations of a noninteracting magnetic monopole. Again, Eq. ()80|) defines a 
dimensionless entity in accord with the fact that the Yang-Mills scale is a parameter 
to be measured and not to be calculated. The right-hand side of Eq. dHU]) measures 
the magnetic flux. If monopoles are point-like, that is, if they are massive, then the 
right-hand side of Eq. (jHUJ vanishes identically due to cancellation of ingoing and 
outgoing fluxes. If monopoles are massless (condensed), that is, if their charges are 
spread over the entire Universe, then this cancellation does not take place, see j33 
for a more detailed investigation. It is clear that definition (|Hn|) relies on definitions 
(jHl) and (j63j) . So when expressed in terms of fundamental caloron and topologically 
trivial fields it looks quite involved. No (nonlocally defined) lattice operator of 




z. m. of n.i. magn. monop. 



(80) 
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this type has ever been constructed. One more point needs to be discussed: The 
phase in Eq. should be a function of Euchdean time r, see Eq. (|H^. From 
the definition of the 't Hooft tensor, Eq. (jU^ . we see that such a time dependence 
manifests itself in terms a r dependent angle in adjoint color space between the fields 
0° and GJJ^. Deep in the electric phase, where monopoles are isolated defects, this 
angle is subject to a global gauge choice for the direction of winding of the field 0'^. 
In the magnetic phase, where monopoles are condensed, the global gauge choice in 
the electric phase is promoted to a local gauge choice for the composite field (p. This 
situation is reminiscent of Kaluza-Klein like gemoetarical compactifications where 
global space-time symmetries along 'extra' dimensions become gauge symmetries 
upon compact ificat ion E2] within a low-energy formulation of the theory. This 
also seems to happens in the low-energy formulation of the Yang-Mills theory being 
in its magnetic phase. 

The fields (pi are energy- and pressure- free if and only if their Euclidean time 
dependence is BPS saturated. Moreover, the ipi must be periodic in time, and their 
gauge invariant modulus must not depend on spacetime. Since BPS monopoles have 
resonant excitations [50J, which are activated by the exchanges of dual gauge bosons, 
one expects the ground-state energy of the system and the tree-level mass spectrum 
of dual gauge-field fluctuations to be T dependent. The local permutation symmetry 
discussed in Secs. l2.6l and l3.1l is respected by the effective potential Vm(v^i, ■ ■ ■ , V^N/2) 
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if it decomposes into a sum over potentials VM{fi)- 

n/2 



The potential Vm is uniquely determined by the above conditions. We have 



VAti^i) = VMi^i)vMi^i) and vm{,Vi) = 'i-Ki/Vi ■ (82) 

In Eq. fl82|) Aa/ denotes a mass scale which is related to the mass scale A^; in the 
electric phase by a matching condition, see Sec. El The effective action for the 
magnetic phase reads 



Sm = J dr I dPx 



(83) 

In Eq. ()83|) G^^y,; denotes the Abelian field strength of the dual field a^i, G^u,i = 
dfj,a^i — dua^i, the covariant derivative is defined as P^,; = dfj, + igaj^i, and g denotes 
the magnetic gauge coupling constant. One remark concerning the normalization of 
the kinetic terms for the fields ifi is in order. The ratio between the gauge kinetic 
term ^ G^u^iG^u^i and the kinetic term | V^ iLpiD^^iipi defines the mass spectrum of 
the gauge-field fluctuations (5a^; in unitary gauge. A redefinition of the factor in 
front of ^^i^iD ^^i^i (and Vm) changes this ratio. At the same time, however, the 
scale Am is changed because the matching condition - equality of the pressure in the 
magnetic and electric phase at the phase boundary - is unchanged. The canonical 
normalization used in Eq. ()83p is thus nothing but a convention for defining the scale 
Km in terms of A^;. 
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The solutions to the BPS equations 

drV^l = VM{<fl) (84) 

read 



A3 

where K{1) is an integer. The condition of MGSB by minimal ground-state energy 
in the electric phase translated into 

K{l) = l. (86) 

Since the l^^ stable monopole species in the electric phase is associated with zeros 
of the SU(2) block in the temporal winding of its condensate ifi is accordingly. 
From Eqs. (|H^ and (jH^ we derive a potential 

V2V-„^MrAi, (87) 

for even N. For N=3 we obtain 

1/2Vm = 'itTAI. (88) 

Again, statistical fluctuations of the fields ipi are negligible and quantum fluctuations 
are absent: 



^^.6,3,3,. (89, 

In Eg. (jHUj) we have defined Am = For N = 2,3 Am is larger than unity 

throughout the magnetic phase, see Sec. 13. 31 Since the fields (pi do not fiuctuate 
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they are a background to the macroscopic gauge-field equations of motion 
There exist pure-gauge solutions to Eq. (jnHjl given as 



(90) 



a^f' = 5,,— T. (91) 

Again, a macroscopic 'holonomy' is created by interacting monopoles which can be 
related to the existence of isolated loops of magnetic fiux: ANO vortices. We have 
i'u,i'^i = 0. As a consequence, the action in Eq. (|H^ reduces to the potential term 
on the ground-state solutions ipi,a^f^. This results in a shift of the ground-state 
energy density and pressure to ±1/2 Vm, respectively. 
Nonperiodic gauge functions 

9i = 27iITt (92) 
transform each pair ipi, a^f^ to unitary gauge 

'^i = \¥^i\, %f^ = 0- (93) 

In analogy to the electric phase one shows that the gauge rotations Qi = e^' leave 
intact the periodicity of the fiuctuations Sa^i- These gauge rotations thus do not 
change the physics upon integrating out the fields 6a^i at one loop. Due to the 
effective theory being Abelian and the monopole condensates (pi inert the one-loop 
calculation is exact. There is, however, a pronounced difference to the electric phase. 
In winding as well as in unitary gauge the Polyakov loop evaluated on each of the 
background fields a^f^ is unity. We conclude that the ground state of the system 
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is unique, much in contast to the electric phase, and thus that the global ^n,o1cc 
symmetry is restored. For a discussion of the full Polyakov loop see Sec. 13.41 Hence 
the magnetic phase confines fundamental test charges at the same time as it allows 
for the propagation of massive, Abelian gauge modes! 

3.3 Gauge-field excitations and 

thermo dynamical self-consistency 

The Abelian Higgs mechanism generates a tree-level mass spectrum for the fluctu- 
ations 6a^i. It is given as 

m = g\'fi\ =aiT . (94) 

Due to the noncondensation of unstable monopoles there are N/2— 1 dual gauge field 
fluctuations 

5a2,, (^ = N/2 + l,---,N) (95) 

which are massless. In analogy to the electric phase we derive an evolution equation 
for temperature as a function of mass from the requirement of thermodynamical 
self-consistency, daP = (for the definition of a in the magnetic phase see Eq. (j98j) ) . 
Notice that the one-loop expression for thermodynamical quantities and the tree- 
level masses of the dual gauge boson fluctuations Sa^f. , (A; = 1, N — 1) are exact 
due to the effective theory being Abelian. We obtain 

96 ""^^ 
" ~ (27r)6N(N+2) ^ ^ ' even) . (96) 
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For N=3 we have 

12 



daKi = --(^^M aDia) . (97) 



In Eqs. dHS} and we have defined: 



_ 1 _ 27rT 

Ci = — , Am = — , 



— \(pi\ = 2ngXM^'^, ai = cia. (98) 



Vr ''"Am 
= -I 

In deriving Eqs. ()9fj|l and (jHTjl we have neglected the 'nonthermal' contribution 
—AVm to the pressure which is very small for sufficiently small N (quantum ffuctu- 
ations are cut off at the compositeness scales \ 

AVm 



<3^jm ~3^X-J. (99) 

1287r2 \Am 1287r2 ^ ^ 



Vm 

The function D{a) is defined in Eq. fjK^ . The evolution equations and (pTTjl 
have fixed points at a = 0, oo which correspond to the highest and lowest attainable 
temperatures in the magnetic phase, Xe,c and Am,c, respectively. 

The Am dependence of the gauge coupling constants g is obtained by inverting 
the solutions to Eqs. and ^7\\ and using the relation between g, Am and a in 
Eq. (j98p afterwards. Results for N=2,3 are shown in Fig.|S| The gauge coupling con- 
stant g increases continuously from g = a.t the electric-magnetic phase boundary 
(Am = ^E,c) until it diverges logarithmically at Am,c- A continuous behavior of the 
magnetic coupling is not in contradiction with magnetic charge conservation since 
no isolated magnetic charges appear in the magnetic phase: magnetic monopoles 
are only present in condensed form, and there are no collective monopole excita- 
tions, see Eq. (jH^ . The continuous rise of g, starting from zero at Am,c, implies a 
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Figure 8: The evolution of the gauge couphng constant g in the magnetic phase 
for N=2 (thick grey hne) and N=3 (thick black line). The gauge coupling constant 
diverges logarithmically, g oc — log(AM — \m,c), at Ajv/,c = 6.41 (N=2) and Xm,c = 
5.82 (N=3). 

continuous behavior of the mass parameter a in Eq. ()98|1 . Since a is the measurable 
order parameter for monopole condensation this situation is reminiscent of a 2'"* 
order phase transition. We compute the critical exponent of this transition for N=2 
in Sec. 13. 51 

Lowering Am towards \m,c the core size of ANO vortices becomes large, the 
monopole condensates are more and more distorted by magnetic flux lines, and thus 
it becomes increasingly irrelevant in what SU(2) embedding a particular monopole 
lives. Formerly unstable monopoles acquire longevity and thus additional monopole 
condensates form at Xm,c- 7^ , (z = N/2 + 1, ■ ■ ■ , N — 1) for Am ~ K,m- As a 
consequence, all dual gauge-boson fluctuations Sa^j, , (fc = 1, ■ ■ ■ , N — 1) are very 
massive close to Am,c and decouple thermodynamically at Am,c- The equation of 
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state at Xm = X^m thus is 



p(Am,c) = -P{Xm,c) 



(100) 



At Xm,c^ where all dual gauge-field fiuctuations are very massive, the continuous 
dual gauge symmetry [/(l)^'^ is broken completely. 



3.4 Polyakov loop in the electric and the magnetic phase 

In this section we show that the Polyakov loop, which is an order parameter for 
the deconfining transition, indeed is finite in the electric and close to zero in the 
magnetic phase. In each phase the Polyakov loop of the full effective theory (now 
we also consider the fluctuations 5ap) formulated in Euclidean spacetime and unitary 
gauge is defined as 

r 2^—1 

P = Z"^ X exp[-^e/] X j VSbp exp -i^ drdh^^- 



exp 



(101) 



where 7 = {e, g} and 6bp = {6ap, Sa^} depending on whether we discuss the electric 
or the magnetic phase. The term exp[— S'cz] refers to the vanishing weight of the 
ground state in the partition function Z in either phase. This weight, however, 
cancels in expectation values. Since the fiuctuations 6bp are periodic in time they 
are decomposed as 



6bp{T,x) 



271171- 



T 



bp,n{x) 



(102) 



n=— 00 

Modes with n 7^ make the Polyakov-loop phase in Eq. (|10ip unity, are action- 
suppressed, and thus they are irrelevant. Zero modes {n = 0) make a contribution 
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Figure 9: The data set used for the fit of the critical exponent v. 

to the Polyakov-loop phase if they are not action-suppressed. This is the case 
if both of the following conditions are met: (i) there is no mass term for these 
fluctuations and (ii) we have dibpfi{x) ~ where di denotes a spatial derivative. In 
the electric phase TLM modes are massless and thus their space-indepent zero-mode 
fluctuations generate the bulk of the (finite) Polyakov loop. For N=2,3 there are 
no massless fluctuations in the magnetic phase and thus condition (i) is violated. 
As a consequence, none of the fluctuations 5bp can contribute to the Polyakov loop 
in a substantial way in the magnetic phase and thus we have P ~ 0. For N>3 
the presence of unstable magnetic monopoles in the magnetic phase prevents some 
of the TLM modes to pick up a mass by the Abelian Higgs mechanism. Thus the 
Polyakov loop should be small but nonvanishing in the magnetic phase. 
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3.5 Critical exponent for the SU(2) electric-magnetic tran- 
sition 

In this section we compute the critical exponent v for the electric-magnetic transition 
for N=2. The obvious order parameter for this transition is the 'photon' mass. 

The data set in Fig.El is generated from an inversion of the numerical solution 
to the evolution equation Eq. The following model is used to fit the data 

a{\M) = Cx \Xm - AA/,e„^r (103) 

where C and u are constants, XM,em denotes the critical temperature Te^c in units of 
and a is the dimensionless 'photon' mass. Recall, that the value Xe^c = 11-65 is 
obtained from the position of the logarithmic pole of the coupling constant e in the 
electric phase. By matching the pressure at the electric-magnetic phase boundary, 
see SecEl this translates into a value XM,em = 7.337. 

To perform the actual fit we have used Mathematica's NonlinearFit function 
which is contained in the statistics package. In Fig. the critical behavior of the 
'photon' mass a is shown. To determine the length of the fitting interval A = 
|AAf,min,flt — ^M,em\ where z/ is least sensitive to changes in AAf,min,flt we numerically 
determine the inflexion point Ainfl^x of the function u = z^(A), see Fig. ^2 We obtain 

Ai„fl,, = 0.29 ± 0.05 . (104) 

We emphasize that this interval is well contained in the fitting interval used in [HH] 
where the critical exponent was determined from the expectation of the dual string 
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Figure 10: The function a{XM) in the vicinity of the electric-magnetic phase tran- 
sition. The region between the dashed vertical lines corresponds to the data set in 
Tab. 1 which generates the least sensitivity of u on the length of the fitting interval 



^ — |AM,min,flt ^ XM,em\- 



tension. Their fit interval < t < 1 with t = 



correponds to an interval 



length A = XM,em ~ 7.34! 

The interval of least senstivity Ai^A^^ = 0.29 ± 0.05 translates into 

z/ = 0.61 + 0.02-0.01. 



(105) 



Alternatively, we can determine z/ by a fit in very small intervals around \M,em- In 
this case we obtain the result expected from a naive mean-field analysis, v —>■ 0.5. 
However, the fitted prefactor C varies considerably for very small intervals around 
^Mfim- It is worth mentioning at this point that the 'would-be' critical exponent 
for N=3 has at least two inflexion points as a function of A. This makes a unique 
determination of the physical value of v impossible and clearly indicates that the 
phase transition is not second order anymore. A small latent heat is associated with 
the electric-magnetic transition for N=3, see Fig.^J which makes it weakly first 
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order. This is seen on the lattice |68j . 

The critical exponent for the 3D Ising model (same universality class as SU(2) 
Yang- Mills [HZIEH]) is z^ising ~ 0.63. As it seems, the effective theories for the electric 
and the magnetic phases have passed an important test! 

4 The center phase 

4.1 Isolated center vortices in the magnetic phase 

Away from the points Xe^c and \m,c there are isolated vortices in the magnetic phase 
which form closed loops due to the conservation of magnetic flux. Along the core 
region of a vortex, where the monopole condensate vanishes, ~ 0, monopoles and 
antimonopoles form a closed chain and move into opposite directions [HSl- Along 
a vortex loop the magnetic flux is 27t k/N , (/c = l,---,N — 1) with respect to the 
dual gauge field aj^j^. This coins the name center vortex loop. The typical action 
S'ano and energy -Eano of a center vortex loop can be estimated since an analytical 
solution to the Abelian Higgs model is known j54j : 

^ANO ~ \ , -^ANO oc - . (106) 

9 9 

As a consequence of Eq. (jl06p center vortex loops do not carry action and are mass- 
less at Am,c where g diverges. They condense to form a new ground state of the 
system: a phase transition takes place. As we will show below this phase transition 
is nonthermal and of the Hagedorn type. 

64 



inflexion 







. 1 



0.2 



0.3 I 
inflex 



. 4 



. 5 



. 6 



A 



Figure 11: Dependence of the critical exponent v on the length of the fitting interval 
A = |AA/,min — Aj\/,c|- The inflexion point of the curve is at AjnA^^ = 0.29 ± 0.05 
corresponding to v = 0.61 + 0.02 — 0.01. 

4.2 Center vortex condensates, macroscopically 

Since center vortex loops are extended, one-dimensional objects the local scalar fields 
$fc(x) , (A; = 1, ■ ■ ■ , N — 1) describing their respective condensates have to be defined 
in a nonlocal way. We formally define the fields in terms of an average over 

the dual Abelian gauge fields A^^^ of the magnetic phase (the part belonging to an 
ANO or center vortex is included in ^^^1) as 



In Eq. fll07j) the integration contour is spatial and circular, its center is at x, and 
circle's diameter is infinite. In absence of a Yang-Mills scale Ac, which is a relevant 
situation when investigating the ground-state structure in the center phase due to 
the missing gauge-mode propagation, only dimensionless quantities like in Eq. ()1()7|1 
can be computed. As we shall see below, the presence of Ac can only be observed in 
the excitation spectrum (selfintersections of center- vortices) . We may always chose 




(107) 
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Figure 12: The phase of a local magnetic center transformation along a circular 
contour parametrized by < ^ < 1 and its decomposition into boxes. 

a parametrization of the integration contour z{^) in Eq. ()l()7j) such that < ^ < 1. 
A (local) magnetic center rotation can be expressed as 



Q{z) = exp 



k=l 



(108) 



where Xk{.z) is either k or zero. The dual gauge field A^j^ transforms under VL{z) as 

K,k^<k-9-'^d,Xk{z). (109) 

According to the definition ()107p a magnetic center rotation VL[z) may locally add 
a magnetic flux quantum [HS] to the flux contained in This is possible since we 
may have VL{z{Q)) ^ VL{z{1)). The action of the local magnetic center rotation VL{z) 
on the complex field $fc(x) is as follows: 

f exp[±2n^] $,(a;), {\xk{z{l)) - x,(^(0))| = k) , 
M^) ^ ■ (110) 

I ^k{x) , otherwise 

Obviously, the action of Q{z) on $fc(a;) may locally change the phase of the field 
and thus the ground state is not invariant under local magnetic center transforma- 
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tions: %^niag 



as a discrete gauge symmetry is spontaneously broken in the center 



phase. 



The local center transformation fl singles out possible 'boundaries of the circle' 
at the positions z{C,n), (^o = 0, ■■■), where it jumps, see Fig.[T21 The size and 
direction of an f2-induced magnetic flux quantum - an observable quantity - should 
not depend on the admissible re-parametrizations^(C), {0 < ( < 1) with^(O) G {^n}- 
For example, a translation 



would have shifted the 'boundary of the circle' at = 0) to = ^i). Thus 
it would have generated a different flux quantum. To avoid such an ambiguity 
we have to impose that a reparametrization of the circle is compensated for by 
an according local permutation of the fields if flux quanta ^ and 

{I k), were generated by Q in the parameterizations ^ and C, respectively, then a 
permutation with <l>; — * needs to be performed after the reparametrization ( was 
implemented. This yields the same flux quantum as in parametrization ^. Now, the 
discontinuous phase change in Eq. piO|) is provided by the dynamics in an effective 
theory and is not externally imposed. As a consequence, the demand for invariance 
of a generated flux quantum under admissible contour-reparametrizations translates 
into an invariance of the effective Lagrangian under local permutations of the set 
■ ■ ■ , <I'n-i(3;)}. The following Lagrangian satisfies this requirement: 



m = c + 6 



(111) 




k=l 



(112) 
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Figure 13: The potential Vc = vc{^)vc{^) corresponding to the definition in 
Eq. pi3|l for N=2,3,4 and Ac = A'(^. |$| is given in units of Ac and Vc in units of 
A^. Notice the minima Vc = at the N"" unit roots. 

It is clear from Eq. ()112|) that the transformation in Eq. ()110|) is only a symmetry of 
the potential term due to the noninvariance of the kinetic terms under the jumps 
in the fields However, because of the conservation of magnetic flux (only closed 
loops of center vortex flux can be generated) one quantum of center flux created 
by a jump at one point in spacetime is compensated by an opposite quantum of 
center flux created by the opposite jump at another point. The spacetime integral 
over Cc, the action, is therefore invariant under local center rotations which induce 
magnetic flux in a closed flux line. Based on the continuum Lagrangian in Eq. pi2|l 
it would be interesting to perform a lattice simulation of the transition to the center 
phase taking any member of the set {|$fe|} as an order parameter. The function 
Vc($) in Eq. ()112p can be written as Vc($) = vc{^)vc{^)- Let us now discuss 
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Figure 14: Numerical solutions to the BPS equation 9^$ = vc{^) for N=2,3,4. 

the properties of the function vc- If a matching to the magnetic phase would take 
place in thermodynamical equilibrium and a classical treatment of the ground-state 
dynamics could be justified at this point then the Euclidean time dependence of the 
(periodic) fields would have to be BPS saturated. Only in this case do the fields 
$fc describe the condensates of zero-energy center- vortex loops 

At first sight a candidate for vc would be fc,triai($fc) = i^%/^k- The potential 
Vc would then not only be ^N,mag but also U{lY'^ symmetric. The latter (global) 
symmetry, however, does not exist in SU(N) Yang-Mills theory and thus fc,triai is 
excluded. 

A function vc{^k)i which is covariant only under transformations and, at 

the same time, allows for periodic and BPS saturated solutions along the Euclidean 

time coordinate r [1^1 , see Fig.lHl is uniquely given as 

^^Recall, that all gauge-boson fluctuations are decoupled in the center phase. As a consequence, 

interaction between center vortices are extremely local. 
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where Ac and A'^- denote mass scales that are a priori independent. The N degen- 
erate minima of the potential Vc{^k) are at 



A^(A'c)H'^". (114) 



At these minima the potential Vc vanishes. Periodic solutions to the BPS equations 

d^^k=vc(^k) (115) 

are parameterized by a winding number G Z in analogy to the situation in the 

magnetic phase. The fields , (/ = 1, ■ ■ ■ , N/2), which are associated with vortices 

formed from the stable dipoles of the electric phase, are winding accordingly. For 

the fields $j , [i = N/2 + 1, ■ ■ ■ ,N — 1), being associated with vortices formed 

from independent monopoles, which are unstable in the electric phase, no winding 

numbers can be derived. This uncertainty in assigning winding numbers for the 

fields $j is refiected in the appearance of an additional mass scale A'^- in Eq. pi3p . 

In the cases N=2,3, however, only stable, independent monopoles exist, and thus 

we can set A^ = A'^^. 

In Fig.HSlthe graphs of the potential Vc for N=2,3,4 and for Ac = A'c are shown. 

Notice the ridges and the valleys of negative and positive tangential curvature, 

respectively. Due to the BPS saturation a solution $fc to Eq. ()115|) is guaranteed 
^•^The latter requirement derives from a consideration of the limit N — > oo where these solutions 

are of physical relevance, see below. 
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to carry no energy. This statement is, however, only then useful for the description 
of the ground state if the classical approximation can be justified. At finite N the 
modulus of BPS saturated solutions is no longer r independent, see Fig.[T31 On the 
one hand, this situation is in contradiction to thermal equilibrium. On the other 
hand, we can not trust the classical approximation since tangential fluctuations Ok-, 
defined as 



can be tachyonic and therefore destabilize the classical solution to Eq. ()115|) . see 
Fig-El In the electric and the magnetic phase tangential fluctuations of the caloron 
and the monopole condensates are would-be Goldstone modes giving rise to longitu- 
dinal polarizations of gauge boson fluctuations. Due to the absence of a continuous 
gauge symmetry in the center phase no gauge-field fluctuations exist which could 
'eat' the tangential fluctuations. How to integrate out tachyonic tangential modes 
analytically is conceptually unclear. The only definite statement we can make at 
present is that they rapidly drive the expectation of the fields towards the min- 
ima of Vc'- relaxes to one of the minima of Vc along an outward directed spiral in 
the complex plane. During this process magnetic flux quanta are locally generated 
by discontinuous phase changes of $fc due the tunneling through the regions where 
the tangential fluctuations are tachyonic. Once the field has settled into one of 
the minima of Vc the situation is classical again. At the minima we have 



$fc = |$fc|exp[i--^] , 




■c 




k 



k 
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Figure 15: The ratio dl^Vc{^k) l\^k? {Vc defined in Eq. (im and Ac 

as a function of |^ for N=2 (thick grey fine), N=3 (thick black fine), N=4 
(dashed fine), and N=10 (thin sohd fine). 

The approximation Ac = A',^ used in Eq. flllTj) is exact for N=2,3. Since radial and 
tangential quantum fiuctuations are heavier than the compositeness scale they 
are absent in the effective theory. As a consequence, the vanishing value of Vc at 
the minima receives no radiative corrections. 

In the limit N — > oo only the pole term of the potential Vc survives for < 
|$™°| and thus we recover the situation of a global U{1)^~^ symmetry. The r 
dependence of the solutions to the BPS equation ()115|) is then a pure phase and no 
tachyonic but only massless tangential fluctuations exist. These fluctuations lead to 
an instantaneous reheating at the phase boundary. The order parameter |$fc| jumps 
from zero to a flnite value across the phase boundary. While the energy density 
of the ground-state jumps to zero at the phase boundary (BPS saturation), the 
ground-state pressure jumps to zero only after = is reached. 

At flnite N a flux quantum and a radial displacement A|$^.| are generated in 
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Figure 16: The creation of an isolated Z2 monopole by self-intersection of a center 
vortex. 

each tunneling process. This reduces the ground-state energy density locally and 
brings the field closer to one of the minima of Vq- The closer $fc to a minimum 
the less likely a tunneling process since the associated energy gain is quickly reduced 
(|5|$^|Vc| decreases!). If the number of tunneling processes taking place until relax- 
ation would be roughly independent of N then for small N many unit-root directions 
will be reached in a relaxation process and thus it is likely that the average value 
of $fc close to the phase boundary is close to zero. For large N only a small an- 
gular sector would be reached by tunneling and thus close to the phase boundary 
the average value of $fc is not close to zero: a strong discontinuity of across the 
transition should be observed. 

Closed center-vortex flux lines can self-intersect and in this way form isolated 
monopoles which reverse the flux, see Fig. for the N=2 case. Let us discuss 
the case N=2 more explicitly. During the phase transition the process leading to 
self intersection can be performed L times on a nonintersecting vortex loop. This 
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generates L isolated Z2 monopoles, each of mass ~ A^, such that the mass 
of the L-monopole state is ~ LKc- Since the number of possible L-monopole 
states roughly grows as L\ (see jnHl for a more precise estimate) we conclude that 
the density of states p{E) for the particle excitations in the center phase of an SU(2) 
theory is over-exponentially growing in energy E: 

E 



p{E) > Ac exp 
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Ac- 

For N>2 the situation is similar. So we conclude that there exists a highest tempera- 
ture Th ~ Ac in the center phase of an SU(N) Yang-Mills theory: a situation which 
was anticipated for any strongly interacting four dimensional field theory a long 
time ago |10j, see also jHO] for a discussion of the Hagedorn transition in an SU(N) 
matrix model. Moreover, we conclude that the center-magnetic phase transition can 
by no means be thermal (the spatial homogeneity of the system is violated during 
the transition) and thus that the thermodynamical pressure may jump across the 
transition. A single and a one-time self-intersecting center-vortex loop for N=2 are 
fermions [59J. This result is crucial for our understanding of the nature of leptons in 
the present Standard Model of Particle Physics. The time scale for the relaxation to 
one of the minima of Vc is roughly given by |($^'''')|^^. A quantitative investigation 
of this reheating process would need methods of nonequilibrium field theory, see for 
example jHI]. In a thermal approach it would be interesting to check in a lattice 
simulation of the effective theory in the center phase Eq. ()112j) whether a density of 
states as in Eq. (jll8|) is indeed seen. 
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The fact that the pressure and the energy density of the ground state are precisely 
vanishing at the minima ^l""" and that there are no radiative corrections to this 
situation clearly is of cosmological relevance. 

5 Scale matching 

The scales A^; and Am, which determine the magnitudes of the adjoint Higgs field 
and the monopole condensates (pi at a given temperature in the electric and the 
magnetic phase, can be related by imposing the condition that the thermodynami- 
cal pressure P be continuous across the thermal electric-magnetic phase transition. 
Disregarding the mismatch in the number of polarizations for some TLM modes in 
the electric phase and some dual gauge bosons in the magnetic phase, we derive 

Ae = {I/A^/'^Am (N even). (119) 

For N=3 we have A^; = (1/2)^/'^ Ajv/. The match between the magnetic and the 
center phase is less determined for the following reasons: (i) For N>3 the winding 
numbers of the fields $j {i = N/2 + 1, ■ ■ ■ , N — 1) close to the center-magnetic phase 
boundary are dynamically generated during the phase transition and so cannot be 
derived from the boundary conditions for the electric phase. For N=2 and N=3 the 
winding numbers of $1 and ($1, $2) are unity, and we can set Ac = A'^; in Eq. ()113p . 
(ii) An analytical description based on the potential Vc of the center-magnetic phase 
transition at finite N breaks down at the phase boundary, see Sec. (Pj). 

In the electric and in the magnetic phase (recall that the latter confines funda- 
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mental test charges) we encounter the thermodynamical analogue to dimensional 
transmutation in perturbation theory: Assuming MGSB, a single, fixed mass scale 
determines the thermodynamics of the SU(N) Yang- Mills theory in these phases. At 
a given temperature this mass scale can be experimentally inferred from the mass 
spectrum of the gauge boson fiuctuations. 



6 Pressure, energy density, and entropy density 



6.1 Numerical results 



In this section we present our numerical results for one-loop temperature evolutions 
of thermodynamical potentials through the electric and magnetic phase. For the 
actual computations we consider N=2,3 only. 
In the electric phase the pressure P divided by is given as 



P_ _ {2Tlf 

T4 ~ 



2A1 



N(N-l) 



\e /N 



^|2(N-l)P(0) + 3 E P{a,)^ + f(t^ + l]^ 



k=l 



(N even) , 
(120) 



where the function P{a) and the dimensionless mass parameters are defined in 
Eqs. and (jlOI), respectively, and MGSB is assumed. In the magnetic phase we 
have 



P 



\4 



9X4 f ^ N/2 

1^ 2(^-l)P(0)+3i:P(a,) 



1=1 



^N(N + 2) 
lb 



(N even) . 
fl211 



The dimensionless mass parameters ai are defined in Eqs. and Eq. 
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Fierure 17: function of temperature for N=2,3. The horizontal hues denote 

the respective asymptotic values. For N=2 we have Xe,c = 11-65 and for N=3 
Xe,c = 8.08. 

For N=3 we have in the electric and the magnetic phase, respectively: 



p 




r2A| 






_(27r)6 


I' 


(2vr)^ 


12XI.J 


2^4 ~ 




_(27r)6 



{4P(0) + 3 (4 P{a) + 2 P(2a)) } + 2Xe 



(N=3) 



rP(a) + A 



M 



(N=3) . 



(122) 



The mass parameters and a; evolve with temperature according to the (inverted) 
solutions to Eqs. (IKT|l . and (FJ7|l . 

Our results for 7^ as a function of temperature in the electric and magnetic 
phase are shown in Fig. El Notice that the pressure is negative in the electric 
phase close to Xe,c and even more so in the magnetic phase where the ground-state 
strongly dominates the thermodynamics. 

In the electric phase the energy density p divided by is given as 



_P_ 

rp4 



At 



2A| 
(2vr) 



N{N-1) 

2(N-l)p(0) + 3 ^ p{ak) 

k=l 



(N even) , 
(123) 
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Fierure 18: function of temperature for N=3 as obtained on a 16^ x 4 lattice 

using the differential method with a universal two-loop perturbative f3 function 
[HSlinSl and using the integral method (solid line) |nZ|- The figure is taken from [H7] . 



where the function p{a) is defined as 



p(a) = / dxx'^ 1^=^= 

"'o exp(va;^ + a^) — 1 



(124) 



In the magnetic phase we have 
P (2vr)^ 



|^j2(^-l)p(0) + 3Ep(.)[ + A|N(N + 2) 



(N even) . 
(125) 



For N=3 we have in the electric and the magnetic phase, respectively: 



P_ 

P_ 

2^4 



(2 



2A| 
(27r)6 

12A^ 



{4p(0) + 3 (4p(a) + 2 p(2a))} + 2\e 



(N=3) , 



M - 



(27r)^ 



p(a) + Aa/ 



(N=3) . 



(126) 



Our results for ^ as a function of temperature in the electric and magnetic phase 
are shown in Fig.^l Slight discontinuities at Xe^c are seen. This is explained by 
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Fierure 19: function of temperature for N=2,3. The horizontal hues denote 

the respective asymptotic values. 

the mismatch in the number of polarizations of fluctuating gauge bosons across the 
electric-magnetic transition - an approximation to scale-matching - and the fact 
that continuity in P does not imply continuity in p. Again, the energy density is 
dominated by the ground-state contribution in the electric phase close to the electric 
magnetic transition and even more so in the magnetic phase. The entropy density 
S is defined as the derivative of the pressure with respect to temperature: 

S.§. (127) 

Using the thermodynamical relation p = — P, we may write 



Our results for S/T'^ are shown in Fig.|2n| The reasons for the slight discontinuities 

at Xe^c are the same as in the case The entropy density S* is a measure for the 
'mobility' of gauge boson excitations. That S is zero at the critical temperature 
Am,c for the center transition is explained by the fact that all dual gauge bosons 



S_ 

T3 
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acquire an infinite mass there and thus no fluctuating degrees of freedom are left in 
the thermodynamical balance. 



6.2 Comparison with the lattice 

An early lattice measurements of the energy density p and the interaction measure 
A = p — 3P in a pure SU(2) gauge theory were reported on in In that work the 
critical temperature Tc for the deconfinement transition was determined from the 
critical behavior of the Polyakov-loop expectation and the peak position of the spe- 
cific heat using a Wilson action. The function A(T) was extracted by multiplying 
the lattice f3 function with the difference of plaquette expectations at finite and zero 
temperature. This assures that A vanishes as T ^ 0. What is subtracted at finite 
T is, however, not the value A(T = 0) since the associated plaquette expectation 
is multiplied with the value of the (3 function at finite T. Apart from this approx- 
imation, the use of a perturbative (3 function was assumed for all temperatures. 
The lattice results for p differ drastically from our results for temperatures close the 
deconfinement, that is, the electric-magnetic transition. We obtain 

— ~ 1.49 (129) 

where Pb = ^T'^ denotes the Stefan-Boltzmann limit (ideal gas of massless gluons 
with two polarizations). On the lattice, this ratio is measured to be smaller than 
unity. At T ~ 5Te,c we obtain 

P 



PSB 



~ 1.34 (130) 
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Figure 20: ^ as a function of temperature for N=2,3. The horizontal hues denote 
the respective asymptotic values. 



while the lattice measures a ratio of about 0.85. Our asymptotic value^'' is 



P 



PSB 



1.33. 



(131) 



Notice the latent heat released at the electric magnetic transition for N=3 fFig.lT^. 

Our result for the pressure P indicates negative values for T close to Te,c (see 
Fig. [T7j) - much in contrast to the positive values obtained in 1^21 • At T ~ ST^ c 
we obtain 

^ ~ 1.30 (132) 

while the lattice measures (P is extracted from the measured values of A and p) a 
ratio of about 0.88. Our asymptotic value for P is 



- 1.32. (133) 

J^SB T^QATe,c 

^^The asymptotic temperature \E,as — 75 is determined by the boundary condition \e{^) = 1000 
for solving the evolution equations H50|) and (|51|l . 
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The asymptotic values of Eqs. ()131|) and ()133|) are very close to the ratio R of the 
number of polarizations for massive TLH modes and massless TLM modes and the 
number of polarizations for massless gluon modes with two polarizations: 

2X3 + 1X2^4^ (134) 

3x23 ^ ^ 

Indeed, at ~ 75 the mass parameter a, defined in Eq. fli^ . is a ~ 27r|^ ~ 0.053 
and therefore the Boltzmann suppression of TLH modes is small. At extremely 
high temperatures a TLH mode 'remembers' its massiveness at low temperatures in 
terms of an extra polarization coming from a tiny mass which, however, still solves 
the infrared problem of naive perturbation theory. 

For a comparison with N=3 lattice data we use the results obtained with a 
Wilson action in [HSj on the lattice of the largest time extension, A^^ = 8. In the 
vicinity of the transition point Te,c the situation for both p and P is similar as for 
N=2: drastic differences between the lattice measurements and our results occur. 
Again, P comes out negative in our calculation, contradicting the positive values 
obtained in [63j. It should be remarked at this point that a lattice simulation of P, 
which did not rely on the integral method (see below) as it was used in [221 ? has 
seen negative pressure for T not far above Tc^e Moreover, the most negative 

value of P/T'^ ~ —0.5 obtained in |66| very close to the phase transition coincides 
with our result at the electric-magnetic transition, see Figs. IT7I and ITHl 
At T = 5 Te^c we have 

~ 1.38 (135) 



PSB 
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while the lattice measures a ratio of about 0.93. Our asymptotic value for p is 



P 



PSB 

For the pressure P we obtain at T = 5 Te^ 



- 1.37. (136) 

T~8.86Te,, 



~ 1.34 (137) 



PsB 

while the lattice measures a ratio of about 0.97. Our asymptotic value for P is 

P , ^ 

- 1.36. (138) 

T~8.86Tb. 



Psi 

Both asymptotic values in Eqs. ()136|) and ()138p are very close to the ratio of polar- 
izations i? = Y = 1.375. 

According to Fig. ()20p the entropy density ^ vanishes at Tm^c- In our approach 
this reflects the fact that at Tm,c all gauge-field are thermodynamically decoupled 
because of their infinite mass. As a consequence, thermodynamics is entirely deter- 
mined by the ground state. This is not observed in the lattice simulations fMj for 
both N=2,3 where a continuous behavior of ^ across the deconfinement transition 
at Te^c was obtained. 

In jnS] a discontinuous behavior of ^ was observed for N=3 using a Wilson 
action and a perturbative beta function. There is an excellent agreement of their 
result with our result for temperatures ranging Te^c down to Tm^-, compare Figs. 1201 
and EU The almost discontinuous behavior in Fig. pUjl is due to the large rise of 
the magnetic gauge coupling with decreasing temperature. According to Fig. (jS)) 
the 'duration' 6T = iE.f±M^ of the magnetic phase is only 6T ~ 0.1. In a lattice 
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Figure 21: ^ as a function of P obtained in SU(3) lattice gauge theory using 
the differential method and a perturbative beta function. The simulations were 
performed on (a) 16^ x 4, (b) 24^ x 4, (c) 16^ x 6 (open circles) and 20^ x 6 (closed 
circles), and (d) 24^ x 6 lattices. Using the 24^ x 6 lattice, the critical value of P is 
between 5.8875 and 5.90. 

simulation the resolution of such a small temperature interval depends very much 
on the choice of the beta function. In a universal perturbative beta function 
was used which may have lead to interprete the behavior of as discontinuous 
in dependence on temperature. Since ^ is a quantity which is much less sensitive 
to the infrared physics than, say, ^ (the direct contribution of the ground state is 
canceled out in ^) the use of (the universal part of) a perturbative beta function 
in the lattice simulations [02] niay be justified. We stress at this point that lattice 
simulations based on lattice sizes of 1-3 times the inverse Yang-Mills scale are not 
capable of being sensitive to the infrared effects of the theory which have correlation 
lengths of the order of the gauge couplings e and g times the inverse Yang-Mills 
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scale at decoupling (Standard Model physics in the electroweak sector suggest that 

~ 0, ~ 10^ n 

^decoup ^decoup ^ * / * 

The alert reader would object that the thermodynamical relation 

dP = SdT, (139) 

which implies that in a homogeneous thermal system the pressure has to be a mono- 
tonic function of temperature, is violated at the center transition (Tc^m), see Fig. lTTl 
where with decreasing temperature the pressure quickly rises from a negative to a 
value close to zero What is the resolution of this puzzle? There are two anwers. 
First, on a microscopic level, the homogeneity of the system at the point where 
center vortices start to condense is badly violated by the generation of (intersect- 
ing) center-vortex loops through discontinuous and local phase changes of the fields 
$fc(x) in Eq. p07|) . The derivation of Eq. ()139p from the partition function, how- 
ever, relies on homogeneity. Second, assuming homgeneity, one can easily convinces 
oneself that the spectral density p{E) in the center phase is exponentially increasing 
with energy E^^: 

p{E) oc exp[;^] . (140) 

H 

Thus the partition function diverges at T = Th ~ Tc^m, the homogeneous system 

would need an infinite amount of energy per volume to increase the temperature 
^^This effect should be the more pronounced the higher N is |58) . 

^^The number of center vortex loops with L intersections increases stronger than factorially with 
L and the energy of a vortex state with L intersections is c>c LAc, for the SU(2) case see [^0] where 
a counting of the vacuum diagrams in a Ac/)* theory is carried out. 
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beyond Th. We conclude that homogeneity is violated at T = Th. We conclude 
that the SU(N) YM dynamics indeed predicts a violation of the thermodynamical 
relation in Eq. ()139|) at T = Th ~ Tc^m 

What are the possible reasons for the qualitative difference between the results 
obtained in jHHl IMI and jHSl ESI? To avoid the use of derivatives of the bare cou- 
pling, which are multiplying the sum of spatial and time plaquette averages in the 
differential method for the computation of the pressure, the integral method was 
introduced in [HZ]. Using a perturbative beta function in the differential method, 
negative values for the pressure for T close to and a rapid approach of the p and 
P to their respective ideal gas limits were obtained in jnEI- The rapid approach to 
the ideal gas limit in p is also obtained in the present approach, compare Eqs. ()130p . 
(jl31|) and Eqs. (|135p . (jl36|) . At the time when the results in jHSl inH] appeared a 
negative pressure was regarded as unphysical and attributed to the use of a pertur- 
bative beta function. The integral method proposed in [HZ] solely generates positive 
pressure. This method assumes the validity of the thermodynamical limit in the 
lattice simulation. Namely, for large spatial volume the thermodynamical relation 

P^T?^ (141) 

valid for a finite volume V is approximated by 

P = (142) 

so that the pressure equals minus the free energy density. In Eqs. (jl41|) and (jl42|) 
^^The author would like to thank D. T. Son for initiatmg this discussion. 
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Z denotes the partition function. Based on Eq. ()142|) the derivative of the pressure 
with respect to the bare couphng /? = ^ can be expressed as an expectation over 
the sum of spatial and time plaquettes without the beta-function prefactor. Thus 
the pressure can be obtained by an integral over (5 up to an unknown integration 
constant. The latter is chosen in such a way that the pressure vanishes at a temper- 
ature well below T^. Instead of only integrating over minus the sum of spatial and 
time plaquette expectations twice the plaquette expectation at T = was added to 
the integrand in jHSHM] to assure that the pressure vanishes at T = 0. We would like 
to stress that this prescription does not follow from relation (jl42|) . Thus it seems 
to be natural that integral and differential methods generate qualitatively different 
results. The results for P(T) obtained by the integral method show a rather large 
dependence on the cutoff and the time extent of the lattice [HSl ■ We believe that 
this reflects the considerable deviation from the assumed thermodynamical limit for 
so-far available lattice sizes. The problem was addressed in where a correction 
factor r was introduced to relate P obtained with the integral method to P obtained 
with the differential method. For a given A'^^ the factor r was determined from the 
pressure ratio at ^ = 0. Subsequently, this value of r was used at finite coupling g 
to extract the spatial anisotropy coefficient (essentially the beta function) by de- 
manding the equality of the pressure obtained with the integral and the differential 
method. In doing so, twice the plaquette expectation at T = was, again, added to 
minus the sum of spatial and time plaquette expectations in the expression for the 
pressure using the differential method. It may be questioned whether a simple mul- 
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tiplicative correction r does correctly account for finite-size and cutoff effects and, 
if yes, whether it is justified to determine r in the hmit of noninteracting gluons^*^. 

To summarize, we observe a quahtative disagreement between the lattice results 
[021 inS] for p and P, using the integral method, and the results of the present 
approach for temperatures close to Te^c- At T ~ ST^ ^ there is better agreement. 
Although for N=3 our negative values for P in the vicinity of Te,c disagree with 
those obtained with the differential method (we obtain a modulus at Te,c which 
coincides with that of the minimal value obtained by using a perturbative beta- 
function differential method jM])- On the other hand, the entropy density obtained 
in jnni with the differential method agrees well with our results. This is expected 
since the entropy is an infrared insensitive quantity. 

7 Conclusions and Outlook 

We have developed a nonperturbative approach to SU(N) Yang-Mills thermodynam- 
ics which is based on the (self-consistent) assumption that the theory 'condenses' 
SU(2) embedded, BPS saturated topological fluctuations of trivial holonomy at an 
asymptotically high temperature. In we have shown for the SU(2) case the re- 
dundancy of this assumption. The concept for the construction of an effective theory 
based on the above assumption is similar to that applied to the construction of the 

macroscopic field theory for superconductivity |22] • We stress that the effects on 
-"^^The Co--values obtained in this way do not coincide with those obtained in |7(J|. 
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nontrivial topology die off in a power-like fashion in the effective theory (as a func- 
tion of temperature). Thus the perturbatively derived suppression of topologically 
nontrivial field configurations does take place in our effective theory as well. 

We have constructed a (uniquely determined) potential for the thermodynamics 
of an energy- and pressure-free adjoint Higgs background which macroscopically 
describes the collective effects due to noninteracting, trivial-holonomy calorons in 
the only deconfining phase of the theory which we call electric phase for obvious 
reasons. As a consequence, the ground state of the system is described by a BPS 
saturated solution to the field equation of the scalar sector and an associated pure- 
gauge configuration solving the macroscopic equation of motion for trivial-topology 
gauge-field fluctuations. The latter has macroscopic, nontrivial holonomy and thus 
describes the presence of isolated magnetic monopoles being generated from decaying 
nontrivial-holonomy calorons as a result of microscopic interactions between trivial- 
holonomy calorons. These interactions are mediated by the trivial-topology sector 
of the theory. The modulus of falls off as oc T~^/^ and the ground-state pressure 
is only linear in T so that nonperturbative effects (apart from extra polarizations 
for excitations) are irrelevant at asymptotically high temperatures. 

Some of the topologically trivial gauge-field fluctuations are massive on tree-level 
due to the adjoint Higgs mechanism, and the underlying SU(N) gauge symmetry is 
spontaneously broken accordingly. An evolution equation, describing the tempera- 
ture dependence of the effective gauge coupling constant e, was obtained from the 
requirement of thermodynamical self-consistency of the effective theory. The two 
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fixed points of tliis evolution were identified. These fixed points predict the existence 
of a highest and a lowest attainable temperature in the electric phase. Based on 
the evolution e(T) a physical argument was given why caloron 'condensation' in a 
grandly unifying theory must take place at a temperature close to the cutoff scale for 
validity of a local, four dimensional field-theory description. We have investigated 
some aspects of the loop expansion of thermodynamical potentials in the effective 
electric theory. Our conclusion is that the present nonperturbative approach resolves 
the infrared problems associated with the usual, perturbative loop expansion. We 
investigate the two- loop corrections to the pressure for N=2 in ^7\. Two- loop con- 
tributions to the pressure are corrections to the one-loop contributions which range 
within the 0.1% level. Thus as far a bulk thermodynamical quantities are concerned 
the gauge-boson fluctuations in the electric fields are practically noninteracting at 
the expense of some of them being thermal quasiparticles. 

The downward temperature evolution of the effective gauge coupling e has an 
attractor and thus the IR-UV decoupling observed in the underlying theory due to 
renormalizability is recovered in the effective theory. The temperature evolution 
e(T) predicts a transition, driven by the condensation of magnetic monopoles, to a 
phase with less gauge symmetry (magnetic phase, confining heavy fundamental test 
charges). This transition is the deconfinement-confinement transition identified in 
lattice simulations. Due to the typical correlation length in the monopole condensate 
being of the order of Cdccoup x A~\ where A denotes the Yang-Mills scale and Cdccoup ^ 
1, it is very hard (in practice impossible) for lattice simulations performed within the 
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magnetic phase to predict thermodynamical quantities that are infrared sensitive. 

The macroscopic ground-state structure of the magnetic phase is determined in 
analogy to the electric phase, and, as a consequence, some of the residual (dual) 
gauge-field fiuctuations acquire mass by the Abelian Higgs mechanism. The evo- 
lution of the magnetic gauge coupling constant g, again being a consequence of 
thermodynamical self-consistency, predicts the existence of a highest and a lowest 
attainable temperature also for the magnetic phase. At the lowest temperature a 
transition to the center phase, where center-vortex loops are condensed into the 
ground state, takes place. Assuming maximal gauge-symmetry breaking in the elec- 
tric phase, aJJ gauge-field modes decouple thermodynamically at the transition point, 
and at this point the equation of state is maximally negative, P — —p. 

The magnetic-center transiton is of the Hagedorn type and thus nonthermal. 
A remarkable feature of the center phase is that ground-state pressure and energy 
density are precisely zero after a period of rapid reheating has taken place. This 
follows from the shape of the effective potential Vc for the vortex-condensate fields 

which implements the spontaneous breakdown of the local center symmetry 

N,mag" 

In the limit N— > oo analytical access to the center-vortex dynamics is granted 
by the thermodynamical and quantum mechanical stability of the classical solutions 
to the BPS equations for the center-vortex fields. For finite N this stability is lost 
due to the presence of tachyonic modes if the center- vortex condensates are away 
from the minima of their potential. Once the minima are reached, no fluctuations 
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in $fc exist for any N > 2, and thus the result of a vanishing ground-state energy 
and pressure, which is based on the classical analysis, is strictly reliable. 

For N=2,3 the present approach predicts a Stefan-Boltzmann like behavior (with 
additional polarizations) of the thermodynamical potentials pressure, energy density, 
and entropy density at temperatures of about lOTjj c- Throughout the magnetic 
phase we predict a negative equation of state which is in contradiction to lattice 
results for N=2,3 using the integral method. For the (infrared insensitive) entropy 
density at N=3 we obtain excellent agreement with lattice data generated with the 
differential method and a perturbative beta function. 

There are many applications of the approach presented in this paper. In [7T] 
we have proposed that a strongly interacting gauge theory underlying the leptonic 
sector of the Standard Model should be based on the following gauge group 

5f/(2)cMB X 5f/(2)e X 5f/(2)^ X SUi2)r ■■■ . (143) 

In addition, mixing angles for the gauge bosons of one factor with those of the other 
factors at temperatures much higher than the Yang-Mills scale of the last factor 
should be supplied. In Eq. (I143j) the Yang-Mills scale of the first factor is roughly 
given (but can be precisely determined |7T]) by the temperature Tcmb of the cosmic 
microwave background ~ 10~^eV. The Yang-Mills scales of the other factors are 
roughly given by the mass of the corresponding charged leptons. While the CMB- 
scale theory is in its magnetic phase very close to the magnetic-electric transition 
(only there is the dual gauge boson, the photon, massless) the other theories are 
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in their center phases generating the stable leptons as sohtons (neutrino ... single 
center-vortex loop, charged lepton ... center-vortex loop with one self-intersection.). 
The two polarization states of these solitons arise as a consequence of the sponta- 
neously broken, local ^2,mag symmetry. The photon, which is the fluctuating degree 
of freedom in the magnetic phase of the CMB-scale theory, would couple to the 
electric charges and the magnetic moments of these leptons because the gauge dy- 
namics subject to Eq. ()143j) was embedded into the gauge dynamics subject to a 
higher gauge group SU(N) at temperatures larger than the mass of the heaviest 
charged lepton. The same reasoning goes through for the coupling of the photon 
to quarks if we allow for SU(3) factors in Eq. (jl43|) . Given the mixing angles it is 
possible to compute the electric charge of each soliton from the plateau value of the 
gauge-coupling evolution in the electric phase in each factor theory. Since the CMB- 
scale theory is in its magnetic phase at the magnetic-electric transition (a very small 
magnetic gauge coupling constant g) the ground state of the universe at present is 
slightly superconducting: a possible explanation for intergalactic magnetic fields. 
The ground-state energy density due to the CMB-scale theory is about 1% of the 
gravitationally observed value [Z1[7T], no contribution. Recall, that no ground-state 
energy density of pressure is generated by SU(N) Yang- Mills theories in their center 
phases. We believe that the missing part can be linked to a CP violating, additional 
term in the Yang- Mills action [7^], see also for another intelligent (albeit incom- 
plete) way of addressing the cosmological constant 'problem'. That the temperature 
of the Universe is stabilized at T = T^mb follows from the behavior of the energy 
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density p at the electric-magnetic transition, see Fig.^J This is an extraordinary 
useful fact since it allows for our mere existence The decoupled bosons of 
SU(2) 

CMB are stable since they cannot decay into the matter that would arise if 

SU(2) 

CMB would go into its center phase (a disaster for entropy generating individ- 
uals). They are an extraordinarily viable candidate for clustering dark matter (the 
stuff responsible for the anomalous rotation curves of galaxies). 

The Zq and the W± bosons of the Standard Model would be interpreted as 
the thermodynamically decoupled dual and TLH gauge-boson fluctuations of the 
SU(2)e factor in Eq. (jl43|) . One would expect to see heavy gauge bosons Z'q and 
W'^, arising from the factor SU{2)^ in Eq. (I143p . at about ^ ~ 200 times the 
mass of the Zq and the W± bosons, respectively. There is no Buctuating Higgs- 
field in this (stepwise) description of electroweak symmetry breaking. Obviously, a 
lot of the extraordinarily precisely checked features of the Standard Model can not 
be derived at the present stage of development, for example the absence of flavor- 
changing neutral currents on tree-level. Moreover, it is questionable that scattering 
processes like e+e^ e+e^, say, at the Zo resonance can be well understood in a 
thermodynamical framework. It is, however, conspicuous that the total cross section 
of this process deviates substantially from the QED prediction for ^/s ~ mg, \T2\ . 
The apparent structurelessness of a charged lepton as measured for momentum 
transfers away from its mass may be understood by the over-exponentially rising 

density of (instable) states in the center phase of the SU(2)e, SU(2)^, ... Yang-Mills 
^^Biology would be unthinkable with a massive photon. 
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theories in Eq. p43|) . 

To make contact with ultra-relativistic heavy ion coUision our approach to pure 
SU(3) Yang-Mills theory would have to be extended to include fundamentally charged 
fermions. The assumption of rapid thermalization, which underlies the (very suc- 
cessful) hydrodynamical approach to the early stages of an ultrarelativistic heavy- 
ion collision 01, would be explained by rigid correlations in the magnetic phase 
(magnetic monopole condensate) or the electric phase (caloron 'condensate') if the 
ground-state structure of a pure SU(3) gauge theory would not significantly be al- 
tered by the presence of quarks. 

To describe thermalized Quantum ChromodjTiamics one would introduce quarks 
as fundamental fermionic fields ipi where i is a flavor index and the color index is 
implicit. Quarks may couple to the caloron 'condensate' (j) in the electric phase via 
Yukawa terms 

KY^^lJiH, (144) 

i 

and to topologically trivial gauge-field fluctuation 5ap via the usual covariant deriva- 
tive. Due to Eq. (|144|) quarks acquire mass dynamically. The ground-state structure 
in this approach would be the same as the one of a pure SU(3) Yang-Mills theory. 
In addition to gauge-field loops there would be quark loops in the expansion of the 
thermodynamical potentials. Thermodynamical self-consistency would imply a sys- 
tem of two coupled evolution equations whose solutions would predict e(T) and k(T) 
and can be used to compute the temperature dependence of the thermodynamical 
potentials. Presumably, the dynamical quark masses would become large close to 
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the transition to the magnetic phase. This would mean that chiral symmetry is 
dynamically broken. As a consequence, quarks would decouple thermodynamically 
at the phase boundary and be replaced by (relatively strongly interacting) chiral 
Goldstone modes in the magnetic phase. The latter could overcompensate the neg- 
ative pressure generated by the ground state of condensed magnetic monopoles. An 
extension of this approach to the case of finite quark chemical potential fig should 
be relatively straight forward. A more fundamental approach, were quarks arise as 
topological solitons in the center phases of various SU(3) Yang-Mills theories, would 
be much more difficult. 

Due to the dominance of the ground state in the magnetic phase a gauge theory 
for cosmological inflation based on SU(N) Yang-Mills thermodynamics would be a 
natural application. This would be a gauge-theory realization of warm inflation 
|76j . Density perturbations generated in the magnetic phase would be dominated 
by thermal fluctuations. Along these lines an attempt to construct a gauge theory 
for warm inflation was made in 

If the entire matter of the Universe would be described in terms of a 'mother' 
SU(N) Yang-Mills theory, which break into SU(K) factors at Mp, then the energy- 
momentum tensor of the ground state would vanish identically for all those 'daugh- 
ter' theories that are in their center phases now. No cosmological constant is gener- 
ated by the latter. 
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